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Executive  Summary 


The  objective  of  this  research  was  to  investigate  a  methodology  for  the  automatic 
development  of  fuzzy  controllers.  Fuzzy  controllers  are  a  class  of  knowledge-based 
controllers  that  use  fuzzy  logic  to  model  control  actions.  They  are  applicable  in  poorly 
understood  or  complex  non-linear  problems  that  defy  control  theory-based  solutions.  A 
number  of  space  applications  are  amenable  to  solution  using  fuzzy  controllers,  including 
battery  charger  control,  fault  diagnosis,  motion  control,  etc.  The  aim  was  to  develop  a 
generalized  domain  independent  system. 

The  two  principal  components  of  a  fuzzy  controller  are  the  fuzzy  rule-base  and  the  fuzzy 
membership  functions.  The  primary  objective  was  to  develop  an  automatic  tuning 
mechanism  for  the  fuzzy  rule-base  (i.e.,  automatically  tune  the  membership  functions) 
using  the  genetic  algorithm  mechanism.  A  secondary  objective  was  to  investigate  the 
applicability  of  genetic  algorithms  to  automatically  tune  and  generate  both  the  fuzzy  rule- 
base  and  membership  functions.  The  genetic  algorithm  is  a  randomized  search  technique 
that  uses  an  evaluation  fimction  to  guide  the  search  process. 

As  an  initial  study  a  fuzzy  controller  package,  the  TIL-Shell  from  Togai  Inffalogic,  was 
studied.  A  typical  control  problem,  the  pendulum  balancing  problem,  was  implemented 
using  the  TIL-Shell.  Issues  relating  to  the  application  of  genetic  algorithms  to  fuzzy 
controllers  were  examined.  A  modified  version  of  the  genetic  algorithm,  the  real  valued 
genetic  algorithm,  was  implemented. 
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In  the  first  .phase  of  this  project,  the  rule-base  (developed  for  the  pendulum  balancing 
problem)  was  tuned  using  the  genetic  algorithm.  The  aim  was  to  tune  the  fuzzy  rule-base 
by  using  the  genetic  algorithm  to  iteratively  improve  the  membership  fimctions.  In  other 
words,  the  genetic  algorithm  searches  through  the  entire  space  of  membership  functions 
attempting  to  minimize  the  error  in  the  output  of  the  fuzzy  controller.  Excellent  tuning 
results  with  high  precision  were  obtained  through  this  process.  A  drawback  of  this 
approach,  however  is  that  it  presupposes  the  existence  of  a  fuzzy  rule-base. 

In  the  second  phase,  the  approach  was  generalized  to  automatically  derive  both  the  fuzzy 
rule-base  and  the  membership  functions  using  the  genetic  algorithm.  The  genetic 
algorithm  searches  through  a  space  of  all  possible  rules  and  membership  functions  to 
derive  a  fuzzy  controller  with  ideal  performance.  As  previously,  the  genetic  search  is 
guided  by  the  minimization  of  the  error  between  the  expected  output  and  the  actual 
output  produced  by  the  fuzzy  controller.  The  full-featured  system  prototype  was 
developed  and  its  performance  verified  with  good  results  on  both  control  and 
classification  problems.  Though  scale-up  issues  were  not  extensively  investigated  in  this 
effort,  this  approach  promises  good  results  for  large,  poorly  understood,  and  complex 
control  applications. 


1.  Fuzzy  Logic  Controllers 

Fuzzy  logic  has  been  widely  applied  in  industrial  control  applications.  This  has  been  the 
result  of  large  complex  system  applications  that  make  precise  measurements  and 
modeling  difficult.  Knowledge-based  controllers  using  linguistic  expressions  and  terms 
are  therefore  a  viable  alternative.  A  qualitative  modeling  technique  that  can  been  model 
the  human  problem  solving  process  is  fuzzy  logic.  The  fundamental  idea  behind  fuzzy 
logic  controllers  is  the  use  of  fuzzy  logic  to  model  and  develop  industrial  control 
applications.  Unlike  traditional  knowledge-based  systems,  fuzzy  control  does  not  use  the 
iterative  execution  of  instructions,  but  rather  the  execution  of  a  number  of  rules  in 
paurallel  -  fuzzy  inference. 

The  following  represents  a  typical  relationship  (rule)  expressed  in  a  fuzzy  system,  a 
“groimd”  truth  (condition),  and  the  inference  (action)  that  may  be  derived  fi-om  these 
relationships. 

Rule  IF  the  distance  between  cars  is  small,  reduce  speed 

Condition  The  distance  between  cars  is  100  ft 


Action  Reduce  Speed  Considerably 

The  action  is  the  result  of  inference.  In  a  fuzzy  controller,  there  would  be  a  number  of 
applicable  rules  that  require  each  inference  rule  to  be  combined  in  some  marmer  and  the 
result  converted  to  a  numerical  value.  This  numerical  value  is  then  used  to  control  the 
system. 
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There  are  four  classical  methods  of  designing  a  fuzzy  controller  [Terano]: 

•  Extraction  of  Domain  Information  from  Experts:  The  experience  of  skilled  operators 
or  the  knowledge  of  control  engineers  is  assimilated  and  brought  together  in  the  form 
of  IF-THEN  rules. 

•  Building  Operator  Models:  In  cases  where  the  skill  of  operators  is  hard  to  quantify, 
their  actions  are  modeled  using  IF-THEN  rules. 

•  Deriving  Rules  Through  Learning:  If  a  plant  model  exists,  fuzzy  rules  can  be  learned 
from  an  environment  in  which  there  are  no  experts. 

•  Development  of  a  Fuzzy  Model  of  the  Plant:  A  fuzzy  model  of  the  system  is  built  and 
fuzzy  control  rules  are  derived  from  the  control  objectives  and  the  plant  model. 

1.1  Fuzzy  Controller  Development 

The  principal  components  of  the  fuzzy  controller  are  the  fuzzy  rule  base,  the  inference 

mechanism,  and  the  input  and  output  interfaces.  The  rule  base  defines  the  actions  of  the 

fiizzy  controller.  There  are  four  steps  in  constructing  a  fuzzy  rule  base: 

•  Identify  and  name  the  input  variables  and  their  ranges. 

•  Identify  and  name  the  output  variables  and  their  ranges. 

•  Define  a  fuzzy  membership  function  for  each  of  the  input  and  output  variables.  The 
membership  functions  assign  a  value  of  inclusion  to  every  possible  input  value.  The 
inclusion  is  with  respect  to  predefined  linguistic  ranges  that  cover  the  entire  or  parts 
of  the  input  and  output  ranges. 

•  Construct  a  rule  base  that  will  govern  the  controllers  operation. 
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•  Determine  how  the  control  actions  will  be  combined  to  form  the  executed  action 


(defuzzification). 

A  number  of  alternate  defuzzification  techniques  have  been  proposed.  In  this  project,  we 
utilized  a  Mamdani-Type  Fuzzy  Controller  [Mamdani]  that  uses  centroid  defuzzification. 
The  centroid  defuzzification  rule  is  given  below; 

Output  =  E  Xi*iJ,A[Xi]/E  ij,aM 

where  X/,  is  the  signal  and  JUAfXi],  is  the  membership  of  Xi,  in  the  linguistic  quantifier, 
A. 

Other  commonly  used  defuzzification  schemes  include  the  Takagi-Sugano  model 
[Takagi].  The  difference  between  the  Sugeno  and  the  Mamdani  model  is  that  the  Sugeno 
model  has  a  mathematical  function  as  the  rule  consequent,  whereas  in  the  Mamdani 
model  each  rule  has  a  fiizzy  consequent. 

The  principal  components  of  the  fuzzy  controller  are  shown  in  Figure  1 . 
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Figure  1 :  Block  Diagram  of  a  Fuzzy  Controller 


The  centroid  defuzzification  algorithm  is  implemented  as  below.  Note  that  a  fuzzy 
controller  simulation  has  to  be  implemented  as  part  of  the  genetic  algorithm  evaluation 


function. 


void  calcMemAreaCent(niember,  lower,  upper.  Area,  Centroid) 

struct  memberObject  *member; 

double  lower; 

double  upper; 

double  *Area; 

double  *Centroid; 

{ 

double  AreaOne; 
double  AreaTwo; 
double  CentroidOne; 
double  CentroidTwo; 

switch(member  ->  type) 

{ 

case  TRIANGULAR: 

♦Area  =  0.5  (member*^  right  Vertex  -  member^  leftVertex)(member->mu); 
CentroidTwo  =  ((niember->rightVertex  -  member">  left  Vertex)  + 
(member^rightVertex  -  member->leftVertex); 

♦Centroid  =  inember->rightVertex  CentroidTwo; 

case  S_FUNCTION: 

AreaOne  ^  0.5  (member^rightVertex  -  member->leftVertex)(member->mu); 
AreaTwo  =  (upper  -  member">centerVertex)^(member->mu); 

♦Area  =  AreaOne  +  AreaTwo; 

CentroidOne  =  member->centerVertex  -  ((member ">centerVertex  - 
member-^  left  Vertex)/3 .0); 

CentroidTwo  =  upper  -  (upper  -  member“^centerVertex)/2.0; 

If  (♦Area  >  0.0) 

{ 

Centroid  =  ((AreaOne^CentroidOne)+(AreaTwo^CentroidTwo))/(^Area); 

}  else 

{ 

♦Centroid  =  0.0; 

} 

case  Z^FUNCTION: 

AreaOne  =  (member->  center  Vertex  -  lower)  ♦(member->mu); 

AreaTwo  =  0.5  ♦(member^  rightVertex  -  member->centerVertex)^(member->mu); 
♦Area  =  AreaOne  +  AreaTwo; 

CentroidTwo  -  member^ centerVertex  -  ((member-^rightVertex  - 
member-^  centerV  ertex)/3 .0); 

CentroidOne  =  lower  +  (member^  centerVertex  -  lower)/2.0; 

If  (♦Area  >  0.0) 

{ 

Centroid  =  ((AreaOne^CentroidOne)+(AreaTwo^CentroidTwo))/(^Area); 

}  else 

{ 

♦Centroid  =  0.0; 

} 

default: 

fprintf(stderr,  “\ncalcMemAreaCent:  Error  in  the  case.  \n”); 
exit(l); 

.} 
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double  CalcOutCentroid(Var) 

/*  this  function  calculated  the  overall  centroid  of  a  variable  */ 

{ 

int  i; 

double  Area,  Centroid,  TotalArea,  TotalCentArea; 

Area  =  0.0; 

Centroid  =  0.0; 

TotalArea  =  0.0; 

TotalCentArea  =  0.0; 

For  (i=0;  i,Var‘^numOfMembers;  i++) 

{ 

calcMemAreaCent(&(Var“>member[i],  Var“>lowerBound,  Var->upperBound,  &Area, 
&Centroid); 

TotalArea  =  TotalArea  +  Area; 

TotalArea  =  TotalCentArea  +  (Area*Centroid); 

} 

if  (TotalArea  >  0.0) 

{ 

return(T  otalCent  Area/T  otalArea); 

} 

else 

{ 

retum(O.O); 

} 


Listing  1 :  Genetic  Algorithm  Evaluation  Function 
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1.2  Developing  Fuzzy  Controllers  with  the  TIL-Shell 

The  TIL-Shell  is  a  software  development  environment  that  provides  a  way  of  describing 
fiizzy  logic  systems,  testing  such  systems  through  simulation,  and  compiling  the  system 
for  further  implementation  on  the  target  processor.  A  fuzzy  system  is  defined  in  the  TIL- 
Shell  using  a  number  of  editors.  V ariables  and  rulebases  are  added  and  connected  in  the 
project  editor.  Membership  functions  are  defined  in  the  variable  or  member  editors,  and 
the  rulebase  is  defined  in  the  rulebase  editor.  The  graphical  project  editor  is  the  top  level 
editor  for  creating  and  manipulating  objects.  It  provides  a  block  diagram  view  of  the 
fuzzy  system. 


1.2.1  Pendulum  Controller 

The  pendulum  controller  is  a  simplified  control  system  for  balancing  a  pendulum.  The 
system  is  developed  using  the  TIL-Shell.  Appendix  A,  shows  the  rule  base  developed 
and  the  membership  functions  defined. 

A  total  of  nine  rules  were  required  for  the  pendulum  controller.  The  rules  are  as  follows: 

IF  error  IS  negative  AND  derror  IS  negative  THEN  current  IS  positive  END 
IF  error  IS  negative  AND  derror  IS  zero  THEN  current  IS  positive  END 
IF  error  IS  negative  AND  derror  IS  positive  THEN  current  IS  zero  END 


IF  error  IS  zero  AND  derror  IS  negative  THEN  current  IS  positive  END 
IF  error  IS  zero  AND  derror  IS  zero  THEN  current  IS  zero  END 
IF  error  IS  zero  AND  derror  IS  positive  THEN  current  IS  negative  END 


IF  error  IS  positive  AND  derror  IS  negative  THEN  current  IS  zero  END 
IF  error  IS  positive  AND  derror  IS  zero  THEN  current  IS  negative  END 
IF  error  IS  positive  AND  derror  IS  positive  THEN  current  IS  negative  END 
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The  definitions  of  the  fuzzy  functions  (in  this  case  positive,  zero,  and  negative,  for  each  of 
the  three  variables,  error,  derror,  and  current)  directly  affect  the  performance  of  the 
controller.  The  tuning  employed  here  consists  of  manually  changing  the  definition  of 
each  fuzzy  function  and  quantifying  the  effect  on  the  controlled  variable.  This  process  is 
laborious  and  infeasible  for  large  controllers. 

2.  Genetic  Algorithms 
2.1  Background 

Genetic  Algorithms  (GAs)  are  a  class  of  iterative,  randomized  search  procedures  capable 
of  adaptive  and  robust  search  over  a  wide  range  of  space  topologies.  They  are  modeled 
after  the  adaptive  emergence  of  biological  species  from  evolutionary  mechanisms 
[Holland].  GAs  have  been  successfully  applied  in  such  diverse  fields  as  image 
processing,  scheduling,  and  engineering  design.  Detailed  introduction  to  GAs  and  their 
applications  can  be  found  in  [Goldberg,  Davis,  Eberhart]. 

The  terminology  used  in  GAs  is  mostly  borrowed  from  the  field  of  genetics.  A  single 
solution  is  called  an  organism  and  the  set  of  solutions  is  termed  a  population.  A 
generation  denotes  a  population  during  an  iteration.  A  string  representation  of  an 
organism  is  called  a  chromosome.  Each  substring  representing  a  feature  of  the  solution 
is  called  a  gene.  The  quality  of  an  organism  is  computed  through  an  evaluation  or 
fitness  function.  Two  mechanisms  of  reproduction  are  used-  1.  crossover,  and  2. 
mutation.  In  crossover,  two  parent  organisms  are  spliced  together  to  create  a  new 


8 


organism.  In  mutation,  arbitrary  bits  in  the  organism  are  flipped.  Mutation  is  generally 
performed  at  a  much  lower  rate  than  crossover. 


The  operational  characteristic  of  GAs  can  be  represented  by  the  following  algorithm: 

•  Create  a  random  population  of  initial  solutions 

•  Compute  the  fitness  of  each  organism 

•  Until  Convergence  Do: 

•  Create  next  generation  by  stochastically  applying 

•  Fitness  Propotional  Selection 

•  Crossover 

•  Mutation 

On  members  of  the  current  population 

•  Compute  fitness  of  members  of  the  new  population 

•  Replace  current  population  by  the  new  population 

•  End  Do 


2.2  Evaluation  of  GA  Packages 

Two  public  domain  GA  packages,  SGA-C  and  Splicer,  were  evaluated  for  robustness, 
ease  of  usage,  and  platform  dependencies.  Splicer  has  a  modular  architecture  that 
includes  a  Genetic  Algorithm  kemal.  Representation  Modules,  and  a  User  Interface 
Libarary.  It  requires  a  UNIX  platform  tvith  X-Windows.  SGA-C  is  a  text-based 
extension  of  the  Simple  GA  program  [Goldberg].  The  primary  platform  is  UNIX. 
Comparative  evaluation  of  the  two  packages  resulted  in  the  selection  of  SGA-C.  Splicer 
had  a  number  of  run-time  problems  associated  with  the  graphical  user  interface.  SGA-C 
had  better  performance  on  sample  problems  and  displayed  greater  robustness. 


Both  SGA-C  and  Splicer  are  based  on  binary  representation  of  the  problem  domain.  The 
mutation  operator  has  significant  effect  on  the  solution.  The  objective  of  mutation  is  to 
achieve  a  local  perturbation  that  moves  the  solution  out  from  a  local  minima.  The  binary 
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nature  of  this  process  creates  large  scale  effects  on  the  solution  rather  than  the  effect  of 
noise.  This  effect  is  known  as  the  binary  hill  problem.  A  solution  to  this  is  to  use  a  real  ' 
valued  GA  (REGA)  where  mutation  is  achieved  by  the  addition  of  noise  to  the  solution. 
Several  previous  studies  have  noted  the  sensitivity  of  membership  functions  and  shapes 
on  the  behavior  of  the  fuzzy  controller.  The  use  of  a  real  valued  GA  could  reduce  the 
sensitivity  of  the  tuning  process.  A  real  valued  GA  was  developed  for  the  UNIX 
platform  (Appendix  B). 

3.  Tuning  Fuzzy  Logic  Controllers  Using  Genetic  Algorithms 

A  number  of  studies  have  shovm  that  Fuzzy  Logic  Controllers  (FLC)  are  sensitive  to  the 
characteristics  of  the  membership  function.  Determining  the  membership  function  is  an  , 
iterative  process  and  requires  engineering  compromises.  Automatic  derivation  of 
membership  functions  can  ease  this  situation  considerably.  In  this  phase  of  the  project, 
we  used  the  genetic  algorithm  paradigm  to  automatically  derive  the  membership 
functions.  Before  the  mechanics  of  this  approach  are  discussed,  the  following  points 
have  to  be  noted: 

•  Fuzzy  sets  of  a  input  or  output  variable  are  ordered  into  regions  of  differing  areas. 
Ideally,  the  area  covered  by  the  fuzzy  sets  decreases  as  the  domain  converges  on  the 
desired  operating  region.  This  has  the  effect  of  reducing  the  number  of  rules 
necessary  to  control  the  system  since  a  single  rule  may  handle  all  the  device  states  in 
the  outlier  regions. 

•  The  degree  of  overlap  at  the  region  of  preferred  performance  provides  fine  tuning  and 
control  by  the  rules.  This  ensures  that  multiple  rules  will  execute  as  the  problem  state 
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moves  t6  the  left  or  right  of  the  operating  region.  For  example,  in  the  pendulum 
balancing  problem,  the  high  degree  of  overlap  ensures  that  tvhen  the  system  is  in  the  .■ 
optimal  state  (the  pendulum  is  balanced),  any  small  changes  are  immediately  detected 

and  handled. 

.  The  number  of  fuzzy  functions  is  always  an  odd  number. 

The  fimdamental  ideal  behind  using  GAs  to  tune  the  FLC  mie  base  is  simple.  Recall  that 
the  GA  attempts  to  converge  to  near  optimal  solutions  using  the  fitness  function.  In  this 
application,  the  fitness  function  is  the  FLC.  The  GA  searches  the  space  of  membership 
functions,  selecting  only  those  memberdiip  functions  that  provide  improved  control.  The 
desimd  input-output  value  pairs  are  provided  to  the  GA.  The  root  mean  square  error 
between  the  output  value  generated  by  the  simulated  FLC  (the  evaluation  function)  and 
the  actual  output  value  is  used  to  drive  the  genetic  algorithm  into  tuning  the  membership 
functions  further.  In  other  words,  we  are  attempting  to  tune  a  black  box  (using  a  GA)  so 
that  it  exhibits  the  desired  characteristics  indicated  by  the  input-output  pans.  For  best 
results  the  input-output  pair  data  should  be  representative  of  the  control  surface. 

Three  classes  of  membership  functions  are  considered  by  the  genetic  algonthms:  the  Z, 
triangular,  and  S  functions.  Any  function  may  be  specified  by  three  vertices:  a  left 
anchor,  a  mid-point,  and  a  right  anchor.  The  GA  searches  domain  of  membership 
functions  (vertices  and  shape)  for  those  that  give  the  lowest  root  mean  square  error 
(indicating  closeness  to  the  supplied  input-output  pairs).  The  rule  base  is  the  set  of  rules 


11 


derived  in  Section  1.2.1.  Figure  2  represents  the  arrangement  of  vertices  in  a 
chromosome. 


IS 

■ 

■ 

■ 

■ 

■1 

■ 

◄-NEGATIVE  ->◄ - zero  - ► 

^ - ERROR - 


LA:  LEFT  ANCHOR 
M  :  MID-POINT 
RA :  RIGHT  ANCHOR 


Figure  2  :  Interpreting  the  Chromosome  -  Fuzzy  Membership  Functions 

3.1  Results 

A  fuzzy  controller  timer  was  developed  using  the  SGA-C  program.  The  pendulum¬ 
balancing  problem  was  used  to  test  the  validity  of  the  controller.  Data  was  gathered  by 
first  defining  the  problem  on  the  TIL-Shell  and  recording  the  input-output  data  pairs  to 
the  pendulum.  Subsequently,  this  data  was  used  as  the  input  to  the  GA-based  tuner. 
There  are  three  components  to  the  GA-based  tuning  process: 

1 .  A  control  rulebase  consisting  of  IF  X  (AND/OR)  Y  THEN  Z. 

2.  A  data  file  consisting  of  input-output  data  pairs. 

3 .  A  GA-based  tuning  mechanism. 
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The  control  rulebase  and  the  fuzzy  membership  functions  together  form  the  evaluation 
function  for  the  GA.  The  control  rule  base  is  the  set  of  rules  derived  in  1.2.1.  Appendix 
C  gives  a  listing  of  the  evaluation  function  (simulated  GA).  Appendix  D  shows  sample 
runs  of  the  tuner.  For  practical  purposes  a  full  run  of  the  GA  has  not  been  included,  but 
the  sample  indicates  the  continuous  improvement  of  the  control  function  (through  the 
modification  of  the  parameters  of  the  membership  functions).  For  instance,  in  Generation 
4,  the  fitness  is  0.430556.  which  by  Generation  199  has  improved  to  0.992738.  The 
maximum  fitness  possible  is  1 .0. 

In  comparison  with  the  manual  system,  the  automated  system  has  two  advantages.  The 
obvious  one  is  the  ease  of  derivation  of  membership  functions.  The  second  is  of 
precision.  It  is  possible  for  the  user  to  specify  the  precision  (number  of  bits)  required  to 
describe  each  input  or  output  membership  function  vertex.  This  gives  higher  precision 
than  the  TIL-Shell,  where  compiler  precision  is  a  limiting  factor.  The  system  that  has 
been  developed  is  very  general.  Any  fuzzy  controller  (irrespective  of  the  domain)  may 
be  tuned.  The  only  requirements  of  the  system  are  a  representative  sample  of  the  input 
and  output  data  pairs  for  the  controller,  and  a  set  of  fuzzy  rules  describing  the  control 
strategy.  This  approach  would  be  most  useful  in  situations  where  the  domain  expert  that 
has  already  generated  the  control  rulebase.  The  evolutionary  approach  eases  the 
difficulty  of  tuning  the  existing  rulebase. 
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4.  Generating  Fuzzy  Controllers  Using  Genetic  Algorithms 

In  the  previous  sections,  we  described  the  development  of  a  fuzzy  controller  tuner  using 
evolutionary  mechanisms.  The  features  of  this  tuner  were  the  following: 

1 .  Evolutionary  learning  of  fuzzy  membership  functions  for  existing  fuzzy  controller 

2.  Tightly  coupled  GA-FCS  mechanism 

3.  Study  of  alternate  GA  mechanisms  to  improve  the  learning  rate. 

The  evolutionary  tuner  was  tested  on  the  pendulum  balancing  problem  and  highly 
accurate  tuning  v/as  achieved.  However,  this  technique  is  limited  to  problems  for  which 
a  fuzzy  rtilebase  already  exists.  It  limits  the  scope  of  applicability  considerably  since  this 
(approach)  presupposes  the  existence  of  such  a  rulebase.  There  are  many  problems  in 
classification,  control,  and  decision  making  for  w'hich  the  structure  of  the  problem  is 
unknown  or  poorly  understood  (i.e.,  the  fuzzy  rules  that  govern  the  control/decision 
process  are  unknown  due  to  the  absence  of  domain  expertise).  This  makes  the 
evolutionary  tuning  process  developed  in  the  first  phase  of  the  project  impracticable.  The 
requirement  here  is  for  a  system  that  can  learn  the  structure  of  the  problem  (i.e.,  the  fuzzy 
rules)  and  achieve  the  tuning  of  these  rules  automatically  (as  in  the  first  phase). 

The  objective  is  to  develop  an  evolutionary  sy.stem  that  could  automatically  learn  the 
structure  of  problems  of  arbitrary  complexity.  The  problem  can  be  subdivided  into  two 
related  sub-problems: 

1 .  Automatic  generation  of  fuzzy  rules 

2.  Tuning  of  the  generated  rules 
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Note,  that  the  prime  objective  is  the  generation  of  articulated  control  systems  (as  opposed 
to  systems  generated  using  neural  networks).  For  the  automated  generation  of  fuzzy- rules 
a  brute  force  teclinique  was  employed.  The  genetic  organism  representation  consists  of 
every  possible  combination  of  inpi-i  and  output  linguistic  quantifiers  (chromosomal 
representation,  one  bit  for  each  rule).  The  assumption  made  was  that  there  are  three 
membership  functions  for  every  input  and  output  variable.  Thus  for  a  system  consisting 
of  m  inputs  and  outputs,  there  are  /'  possible  rules.  The  rule  generation  part  of  the 

organism  will  thus  have  4"'  bits. 

Figure  3  is  a  representation  of  the  GA  organism.  It  is  a  generalization  of  Figure  2,  with 
additional  genetic  material  added  for  the  representation  ot  rules. 


FUNCTIONS 

LA : LEFT  ANCHOR 
M  :  MID-POINT 
RA :  RIGHT  ANCHOR 

Figure  3:  Interpreting  the  Chromosome  -  Fuzzy  Rulebase  and  Membership  Functions 

The  GA  mechanism  is  used  to  search  through  the  space  of  all  possible  rules  and  valid 
membership  functions  to  obtain  the  best  control  system.  An  important  objective  is  rule 
minimization.  Solutions  with  minimal  number  of  rules  are  preferred  to  controllers  with  a 


larger  number.  A  rule  generator  was  designed  and  coded.  The  rules  are  indexed  into  a  2- 
dimensional  linked  list.  The  rules  have  been  integrated  into  the  genetic  algorithm 
environment. 


5.  System  Verification 

The  system  was  exhaustively  tested  on  the  pendulum  balancing  problem  and  a  tank 
recognition  problem.  The  pendulum  balancing  problem  is  a  two  input  -  one  output 
problem.  A  complex  non-linear  relationship  exists  between  the  input  and  output 
variables. 


5.1  Results 

Pendulum  Balancing  Problem  :  The  control  problem  is  a  two  input  (Error  and  Derror)  - 
one  output  (Current)  problem.  The  system  generates  the  set  of  all  possible  rules  that  can 
describe  the  input-output  relationship.  This  rule  space  (generated  automatically  by  the 
system  developed)  is  shown  in  Appendix  E.  There  are  three  membership  functions 
defined  per  variable  -  negative,  zero,  and  positive.  The  rule  space  consists  of  45  rules 
and  9  membership  functions.  The  membership  functions  may  be  S-functions,  Z- 
functions,  or  triangular. 

The  results  of  the  evolutionary  generation  of  the  fuzzy  control  system  are  detailed  in 
Appendix  F.  The  system  has  reduced  the  rule  set  to  six  rules.  The  definition  of  the  fuzzy 
memberships  are  also  given  in  Appendix  F.  The  performance  results  are  very  good  with 
a  Mean  Square  Error  of  0.039. 
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T^nV  Recognition  Problem:  The  problem  considered  is  the  recognition  of  tanks  in  a 
battlefield  based  on  the  audio  sensor  information.  The  data  consists  of  three  inputs  and 
one  output.  A  complex  non-linear  relationship  exists  between  the  input  and  output 
variables;  however,  it  is  not  known  what  the  relationship  is.  A  set  of  apriori  classified 
data  is  used  to  train  and  subsequently  test  the  system.  At  the  start,  unlike  the  pendulum 
balancing  problem  (which  was  manually  coded  and  tested),  it  was  not  known  how  many 
rules  would  be  required  to  describe  the  system.  The  system  generates  the  set  of  all 
possible  rules  that  can  describe  the  input-output  relationship.  This  rule  space  (generated 
automatically  by  the  system  developed)  is  shown  in  Appendix  F.  Note  that  there  are 
three  membership  functions  defined  per  variable  -  negative,  zero,  and  positive.  The  rule 
space  consists  of  189  rules  and  nine  membership  functions.  The  membership  functions 
may  be  S-functions,  Z-functions,  or  triangular.  The  results  of  the  evolutionary  fuzzy 
control  system  are  detailed  in  Appendix  F.  Note  that  the  system  has  reduced  the  rule  set 
to  seven  rules.  The  performance  results  are  reasonably  good  with  a  Mean  Square  Error  of 
0.414  and  a  Floor-Ceiling  Mean  Square  Error  of  0.498.  A  recognition  rate  of  76%  is 
achieved  on  the  test  data.  The  larger  error  is  due  to  the  discretization  of  the  output  of  the 
fuzzy  controller  required  by  the  classification  problem.  Another  issue  is  that  the  full  tank 
recognition  problem  was  not  utilized  (i.e.,  the  entire  input  parameter  set).  This  approach 
attempted  the  classification  process  using  the  three  principal  components  of  the  problem. 
Further  tuning  and  the  use  of  the  entire  parameter  set  can  improve  the  recognition  rate. 
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A  serious  drawback  with  this  approach  is  the  combinatorial  explosion  (of  possible  rules) 
when  the  control  problem  has  a  large  number  of  inputs  and  outputs.  While  the 
combinatorics  is  not  a  problem  for  the  genetic  algorithm,  adequate  representation  of  the 
rule  set  is  problematic.  In  this  implementation,  rules  are  recursively  generated  for  the 
complete  set  of  input  and  output  variables.  Even  small  sets  of  these  variables  can  result 
in  the  computational  overhead  of  very  large  organisms. 

6.  Conclusions  and  Future  Work 

Theoretical  issues  regarding  the  application  of  genetic  algorithms  to  the  development  of 
fuzzy  logic  controllers  were  studied  in  this  research  effort.  This  study  has  resulted  in  the 
development  of  a  significant  research  prototype.  A  fuzzy  rulebase  was  developed  using  a 
commercial  fuzzy  controller  shell  (the  TIL-Shell)  and  tuned  manually  and  automatically. 
The  results  of  the  automatic  tuner  are  good,  and  deliver  much  higher  precision  than  the 
manual  tuner.  There  is  also  a  significant  reduction  in  the  effort  required  in  producing 
robust  controllers.  However,  this  approach  presupposes  the  existence  of  a  fuzzy 
controller  rulebase.  This  is  always  not  the  case.  To  further  generalize  the  approach,  the 
prototype  evolutionary  fuzzy  control  system  was  enhanced  (though  not  part  of  the  initial 
project  plan)  to  generate  both  fuzzy  rules  and  tune  the  membership  functions 
simultaneously. 

In  general  terms,  this  approach  is  applicable  across  a  wide  spectrum  of  practical 
problems.  The  objective  of  this  research  was  to  investigate  space  applications  for  this 
technique.  Two  space  applications  suitable  to  the  sponsor  were  considered.  One  of  these 
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applications  involved  the  control  of  smart  structures.  Data  was  not  immediately  available 
for  this  application.  In  the  other  application  Lockheed  Martin  supplied  controller 
simulation  data  for  control  of  the  null  position  of  a  solar  drive  encoder  and  to  set  the 
encoder  initialization  bits.  The  results  of  the  evolutionary  fuzzy  controller  were 
inconsistent  with  large  mean  square  errors.  Table  1  shows  a  portion  of  the  Lockheed 
Martin  data. 


+Y_Plane_Error 

0.1758 

-1.0 

-4.771 

1 

0.1758 

0.0 

-4.771 

2 

0.1758 

1.0 

-4.771 

3 

0.1758 

2.0 

-4.771 

4 

0.1758 

3.0 

-4.771 

Table  1 :  Solar  Drive  Control  Data 

The  +Y_Step  is  the  output  variable  and  the  others  the  input  variables.  Two  points  about 
this  data  are  noteworthy.  First,  Time  is  a  monotonically  increasing  variable  that 
determines  +Y_Step.  The  evolutionary  GA  controller  does  not  handle  monotonically 
increasing  or  decreasing  variables  well.  Second,  the  output  is  dependant  on  Time  alone. 
The  unsatisfactory  results  produced  by  the  GA-based  controller  are  related  to  the  nature 
of  this  data.  This  application  does  not  lend  itself  to  a  knowledge-based  controller 
solution.  However  this  is  a  feature  of  the  specific  problem  domain  rather  than  a 
drawback  of  the  technique.  Further  interaction  vsith  the  Phillips  Laboratory  would  reveal 
a  wider  class  of  problems  that  can  be  easily  solved  by  this  technique. 
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Currently,  the  prototype  is  capable  of  generating  valid  fuzzy  controllers  in  cases  where 
the  number  of  inputs  and  outputs  is  small.  In  general  good  results  were  obtained  on  both 
classification  and  control  problems.  A  problem  associated  with  the  GA  approach  is  that 
the  space  of  valid  controllers  is  very  small,  and  typical  GA  operations  of  mutation  and 
crossover  can  quickly  move  a  controller  from  a  valid  to  an  invalid  region.  This  does  not 
give  sufficient  time  for  the  GA  to  iteratively  improve  the  solution.  A  solution  could  be 
the  development  of  semantically  driven  mutation  and  crossover  operations  that  narrow 
the  search  space  rapidly. 

As  mentioned  previously  representation  is  a  significant  issue  when  dealing  with  large 
numbers  of  inputs  and  outputs.  We  have  done  preliminary  investigation  on  an  alternate 
methodology  that  overcomes  the  combinatorics  in  the  representation  -  the  Combs 
technique,  with  good  results.  In  the  Combs  technique  every  input  is  mapped  to  every 
output.  The  resulting  system  has  linear  complexity.  Tuning  such  a  system  is  more 
difficult,  but  the  GA-based  tuner  would  ease  this  situation.  This  is  a  fruitful  area  for  the 
investigation  into  complex  fuzzy  system  applications. 
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Appendix  A 

Pendulum  Controller  Developed  Using  the  TIL-Shell 

PROJECT  Project  1 
OPTIONS 

ICONCOLOR=12632256 

MODE=’'NORMAL" 

CHANGEID=2690521 803 
END 


/*  The  membership  functions  were  defined  using  the  membership  editor  */ 
/*  Membership  functions  for  Error  */ 

VAR  Error 
OPTIONS 
ICONPOS=0.5,0.5 
FULLGRAPHICS="ON" 

GRIDSHOW="OFF" 

GRIDSNAP="OFF" 

GRIDSPACE=0.2.0.2 

NUMBER=3 

SCALE=50 

TOUCHED='’ON" 

END 

TYPE  float 
MIN-1 
MAX  1 

MEMBER  N 
OPTIONS 

ICONCOLOR=1671 1680 
END 

POINTS  -1,1  -0.58024691358,0.989583333333  0.037037037037,0 
END 

MEMBER  Z 
OPTIONS 

ICONCOLOR=65407 

END 

POINTS  -0.446913580247,0  0.140740740741,1  0.5,0 
END 

MEMBER  P 
OPTIONS 
ICONCOLOR=255 
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END 

POINTS  0.0172839506173,0  0.264197530864,1  1,1 
END 
END 

/*  Membership  functions  for  dError  */ 

VAR  dError 
OPTIONS 
ICONPOS=0.5,2.5 
GRIDSHOW="OFF'' 

GRIDSNAP="OFF" 

GRIDSPACE=0.2,0.2 

NUMBER=3 

SCALE=50 

TOUCHED="ON" 

END 

TYPE  float 

MIN-1 

MAXI 

MEMBER  N 
OPTIONS 

ICONCOLOR=1671 1680 
END 

POINTS  -1,1  -0.257731958763,1  0.0103092783505,0 
END 

MEMBER  Z 
OPTIONS 

ICONCOLOR=65407 

END 

POINTS  -0.365979381443,0  -0.20618556701,0.987341772152  0.5,0 
END 

MEMBER  P 
OPTIONS 
ICONCOLOR=255 
END 

POINTS  -0.0257731958763,0  0.515463917526,1  1,1 
END 
END 

/*  Membership  functions  for  Current  •/ 

VAR  Current 
OPTIONS 
ICONPOS=3,1.5 


23 


GRIDSHOW=''OFF" 

GRIDSNAP="OFF" 

GRIDSPACE=0.2,0.2 

NUMBER=3 

SCALE=50 

TOUCHED="ON'' 

END 

TYPE  float 
MIN-1 
MAX  1 

MEMBER  N 
OPTIONS 

ICONCOLOR=16711680 

END 

POINTS -1,1  -0.551546391753,1  0.025773 ’.958763,0 
END 

MEMBER  Z 
OPTIONS 

ICONCOLOR=65407 

END 

POINTS  -0.278350515464,0  -0.134020618557,0.987341772152  0.417525773196,0 
END 

MEMBER  P 
OPTIONS 
ICONCOLOR=255 
END 

POINTS  0.0257731958763,0  0.5,1  1,1 
END 
END 


RULEBASE  Pend_Rules 
OPTIONS 
ICONPOS=1.5,1.5 
END 

RULE  Rulel 

IF  (Error  IS  N)  AND  (dError  IS  N)  THEN 
Current  =  P 

END 

RULE  Rule2 
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IF  (Error  IS  Z)  AND  (dError  IS  N)  THEN 
Current  =  P 

END 

RULE  Rules 

IF  (Error  IS  P)  AND  (dError  IS  N)  THEN 
Current  =  Z 

END 

RULE  Rule4 

IF  (Error  IS  N)  AND  (dError  IS  Z)  THEN 
Current  =  P 

END 

RULE  Rules 

IF  (Error  IS  Z)  AND  (dError  IS  Z)  THEN 
Current  =  Z 

END 

RULERule6 

IF  (Error  IS  P)  AND  (dError  IS  Z)  THEN 
Current  =  N 

END 

RULE  Rule? 

IF  (Error  IS  P)  AND  (dError  IS  P)  THEN 
Current  =  N 

END 

RULE  Rules 

IF  (Error  IS  Z)  AND  (dError  IS  P)  THEN 
Current  =  N 

END 

RULERule9 

IF  (Error  IS  N)  AND  (dError  IS  P)  THEN 
Current  =  Z 

END 

END 

DEBUG  Debugl 
END 

CONNECT 
FROM  Error 
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TO  Pend  Rules 
END 

CONNECT 
FROM  dError 
TO  Pend_Rules 
END 

CONNECT 

FROM  Pend  Rules 
TO  Current 
END 
END 
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APPENDIX  B 

REAL  VALUED  GENETIC  ALGORITH  (ReGA) 


/•  Real  Valued  Genetic  Algorithm 
/*  Written  by  :  Radhakrishnan  Srikanth,  Roy  George 
/*  Department  of  Computer  Science 
/♦  Cl^  Atlanta  University 


/******••••♦••*••****•*•*♦•*♦*****♦***•♦*♦****••**♦•********•*****/ 


*/ 

•/ 

*/ 

*/ 


#include  <stdlib.h> 

^include  <stdio.h> 

#include  <math.b> 

^include  <time.h> 

#defineeq  = 

^define  MAXINT  777777 

typedef  float  GeneType; 

typedef  struct  pop  {  I*  define  item  to  hold  the  value  and  the  fitness  *! 

GaieType  ‘chromosome; 
float  fitness; 

}PopType; 

void  InhializeO;  t*  Initializes  the  members  of  the  population  *f 

void  SrsSelectO;  /*  Selects  members  of  the  Next  Generation  */ 

void  CrossOverO;  /•  Selects  and  CrossesOver  Chromosomes  with  a  given  probability  */ 

void  CrossOverChromosomeO;  /*  Crossover  generates  two  new  children  from  given  Parents  */ 

void  MutateO;  /*  Mutates  Chromosomes  in  the  population  •/ 

void  MutateChromosomeO;  /*  Mutates  parts  of  the  chromosome  */ 

void  CalculateFitnessO;  /*  Calculates  fitness  of  the  chromsome  •/ 

void  PrintReportO;  /*  Prints  Report  •/ 

void  PrintChiomosomeO;  /*  Prints  induvidual  Chrcsnosome  */ 

void  CorssOverO;  /*  Performs  Crossover  */ 

void  CopyChromosomeO;  /*  Copies  a  give  chromosome  */ 

float  ObjectiveFunctionO;  /*  Calculates  the  fitness  of  a  chromosome*/ 

void  SRSeiectQ;  /*  Selects  the  new  p(q>ulation  */ 

void  StatisticsO;  /*  Keeps  foe  Statistics  for  each  generation  */ 

void  InitializeDataO;  /*  User  written  routine  which  reads  data*/ 

void  TestFuncticMiOi  Prints  foe  ou^ut  */ 

int  PopulatitxiSize;  /*  Population  Size  */ 

int  NumberOfParameters;  /*  Number  of  Parameters  per  ChromoscHne  */ 
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iot  seed;  . 
random  nximbers*/ 


/*  The  random  seed  that  is  used  for  generating 


float  CrossOverProbability;  /*  Stores  the  CrossOver  Probility  ♦/ 
float  MutateProbability;  /•  Stores  the  Mutation  Probility  */ 

float  range=I ;  /*  Upper  bound  of  range  of  the  mutation  amount  ♦/ 

int  CrossOverStat=0;  /*  Number  of  Crossovers  */ 

int  MutationStat=fl;  /*  Number  of  Mutations  */ 

int  GenerationNumber;  /♦  Keeps  track  of  the  generation  */ 

FILE  '"graph;  /*  Keeps  the  statistics  *! 

PopType  BestEver;  /*  Stores  the  Best  ever  Chromosome  •/ 


FILE  “"Data; 
exan:q)le  */ 


/*  Data  file  for  the  neural  network  training 


int  NumExamples;  /*  Number  of  examples  in  the  file  *! 

int  Numinputs;  /•  Number  of  input  units  in  neurlal  network*/ 

int  NumOu^uts;  /*  Number  of  output  units  in  neurlal  networic*/ 

int  NumHidden;  /*  Number  of  hidden  units  in  neurlal  network*/ 

float  InputPattem[2][4];  /*  Stores  training  pattern  */ 
float  Tai;get[4];  /*  Stores  target  pattern  */ 


void  mainO 

{ 

int  i;  /*  Counters  */ 

*Populationl;  /*  Contains  current  population  of  genes  */ 

PopType  *Population2;  /*  Contains  next  population  of  genes  */ 
int  MaxGenerations;  /*  Maximum  number  of  genertions  */ 

/*  Prompt  and  Read  number  of  parameters  and  the  size  of  the  population  */ 

printf("Enter  the  population  size\n"); 
scanft"%d",&PopulationSi2e); 

printf("Enter  the  number  of  parameters  in  each  Chromosome\n"); 
scanf("%d",&NumberO£Parameters); 

/*  allocate  space  for  the  members  of  the  population  */ 

Population!  =  (PopType  *)  malloc(PopulationSize*sizeof{PopType)); 

Population2  =  (PopType  *)  maIloc(PopulationSize*sizeof(PopType)); 

/*  allocate  space  for  each  chromozome  */ 

for(i=0;i<PopulationSize;i-H-){ 

Population! [i].chromosome  =  (GeneType  *)  malloc(NumberC>fParameters  * 

sizeof(GeneType)); 

Population2[i].chromosome  =  (GeneType  *)  malloc(NumberOfParameters  * 

sizeof(GeneType)); 

} 
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BestEver.chromosome  =  (GeneType  *)  inalloc(NumberO£Paranieters  *  sizeof(GeneType)); 
BestEver.fitness  =  (£loat)(-MAXINT); 

Initialize(Population  1 ) ; 

printfC'Enter  Cross  Over  Probability  (a  number  between  0-l)\n"); 

scanf("%f',&CrossOverProbability); 

printfC'Enter  Mutation  Probability  (a  number  between  0-l)\n"); 

scanf("%f',&MutateProbability); 

printfC'Enter  number  of  generations\n"); 

scanf("%d",&MaxGenerations); 

/•  Open  the  ou^ut  file  */ 

ifi[(grapb  =  fopen("output'',"wt-"))  =  NULL){ 
printfC'Unable  to  open  ou^ut  file\n"); 
exit(0); 

} 

InitializeDataO;  /*Read  Training  Examples  for  the  Neural  Net  */ 


for(GenerationNumber=0;GenerationNumber<MaxGenerations;GenerationNumberH-){ 
CalculateFitnessfPopulation  1 ); 

SRSelectfPopulation  1  ,Population2) ; 
CrossOver(Population2,Population  1 ); 

MutatefPopulation  1 ); 

Statistics(Population  1  ,&BestEver); 


} 

PrintReportfPopulation  1 ); 
fi^Population  1 ); 
fi'ee(Population2); 
fclose(graph); 


} 


void  Initialize(PopType  ’Population) 

{ 

intij; 

printfC'Enter  any  integer  for  a  seed\n'’); 

scanfC'%d",&seed); 

srand(seed); 


15.12); 


for(i=0;i<PopulationSi2e;i-H-)  { 

fijr0=0j<'JumberO£Parametersj-H-) 

Population[i].chromos(Hne[jl  =  (GaieType)(((randO%100)/3.3)  - 


29 


} 


} 


void  CalculateFitaess(PopType  ‘Population) 

{ 

int  i; 

for(i=0;i<PopulationSize;i-H-)  { 

Population[i]. fitness  =  ObjectiveFunction(Population[i].chromosome); 

} 

) 


void  Statistics(PopType  ‘Population,  PopType  ‘BestEver) 

{ 

inti; 

float  Best; 
float  Worst; 
float  Average; 

Best  =  (float)  (-MAXINT); 

Worst  =  (float)  (MAXIMT); 

Average  =  (float)  0; 

for(i=0;i<PopulationSize;i-H-){ 

Population[i]  .fitness  =  ObjectiveFunction(Population[i].chroniosome); 
/‘  ifl[Population[i].fitness  >  BestEver->fitness) 

CopyChiomosome(PopuIation[i],BestEver);  *! 
if(Population[i]  .fitness  >  Best) 

Best  =  Population^]  .fitness; 
if(Popuiation[i]. fitness  <  Worst) 

Worst  =  Population[i]  .fitness; 

Average  +=  Popu)ation[i]  .fitness; 


} 

Average  =  Average/PopulationSize; 

^rintf(graph,"%d  %f  %f\n",GenerationNumber,Best,Worst); 

} 


/‘  Prints  report  ‘/ 

void  PrintReport(PopType  ‘Population) 
intij; 

for(i=0;i<PopulationSize;i++){ 
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printf("\nChromosome  Fitness  %f\n",PopuIation[i]. fitness); 
forO=Oj<NumberOfParametersJ-H-) 

prmtf("Population[%2d].chromosome[%2d] 

%f\n",ij,Population[i].chromosome|j]); 

printfl("\n''); 

} 

printfi^’Total  number  of  Crossovers  %d\n",CrossOverStat); 
printfC’Total  number  of  Mutations  %d\n”,MutationStat); 
printfC'Best  Ever  fitness  is  %f\n",BestEver.fitness); 
printf(”Best  Ever  Chromosome  is:  \n"); 

/*  fbrO=Oj<NumberO£ParametersJ-H-) 

printf("Population[%2d]  .chromosome[%2d] 
%f\n",ij,Population[i].chromosome|j]); 

•/ 


} 


TestFunction(BestEver.chromosome); 


void  PrmtChromosome(GeneType  *chromosome) 

{ 

intj; 

for(j=Oj<NumberOfParametersJ-H-) 

printf(''chromosome[%2d]  =  %f\n"j,chromosome|J]); 


} 


/*  InitializeData  function  written  by  user  to  read  appropriate  Global  Data  *! 

void  InidalizeDataO 

{ 

intij; 


if((Data  =  fopen("input.dat","r"))  =  NULL){ 

printfC’Error  opening  file  input\n"); 
exit(-l); 

} 


&canf(Data,"%d  %d  %d",&NumInputs,&NumHidden,&NumOutputs); 

&canf(Data,"%d",&NumExamples); 

for(i=0;i<NumExamples;i-H-)  { 

for(j=0  J<NumInputs;j-H-)  { 

&canf(Data,"%d",&lnputPattem{j][i]); 

} 

fecan£(Data,"%d",&Target[i]); 


} 
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fclose(Data); 

} 


/*  Objective  function  typically  written  by  the  user  */ 

float  ObjectiveFunction(GeneType  *chromosome) 

{ 

int  y,k,l; 
float  net; 
float  out; 
float  hout[2];  . 

float  fitness;  /*  fitness  of  the  chromosome  */ 
fitness  =  (float)0.0; 

/*  for  (i=0;i<NumberO£Parameters;i-H-) 

fitness  +=  chromosome^];*/ 

for(i=0;i<NumExamples;i-H-)  { 
k=0; 

for(l=0;l<NumHidden;l-H-)  { 
net  =  (float)0.0; 
for(j=0  J<NumInputs  J-H-)  { 

net  +=  InputPattem[i][j]*chromosome[k]; 

k-H-; 

} 

net  =  net-chromosome[k]; 
k++; 

houtpj  =  (float)net/(l+febs(net)); 

} 

for(l=0;l<NuniHidden;l-H-)  { 

out  +=  hout[l]*chromosome[k]; 
k++; 

} 

out  =  out  -  chromos(Mne[k]; 

k-H-; 

out  =  (float)out/(l+fiibs(out)); 
if(fiibs(out-Target[i])  <  0.5){ 
fitness  +=  (float)  1.0; 

} 

} 

ietum(fitness); 


} 
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void  TestFunction(GeneType  •chromosome) 

{ 


int 

float  net; 
float  out; 
float  bout[2]; 

for(i=0;i<NumExamples;i-H-)  { 
k=0; 

forO=0;l<NumHidden;l-H-)  { 

net  =  (float)0.0; 

for(j=0  j<NumInputs  J-H-)  { 

net  +=  Ii^utPattem[iJ[j]*chromosome[k]; 
k++; 

} 

net  =  net-chromosome[k]; 
k-H-; 

houtp]  =  (float)net/(l+febs(net)); 

} 

for(l=0;l<NumHidden;l-H-){ 

out  +=  hout[l]*chromosome[k]; 
k-H-; 

} 

out  =  out  -  chromosome[k]; 
k-H-; 

out  =  (float)out/(l-<-febs(out)); 
printfC’output  for  pattern  %d  is  %f\n",i,out); 


} 

/•  CrossOverChromosome,  is  a  disruptive  operator  will  take  two 
induviduals  in  the  populaticm  and  produce  two  children  •/ 


void  CrossOverChromosome(GaieType  •parent  l,GeneType 

•child  l.GeneType  •child2) 

{ 

int  CrossOverPoint;  /•  Randomly  Chosen  Cross  Over  Point  •/ 
int  i;  /•  Counters  •/ 


*parent2,GeneType 
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void  TestFunction(GeneType  •chromosome) 

{ 


int  i 

float  net; 
float  out; 
float  bout[2]; 

for(i=0;i<NumExamples;i++)  { 
k=0; 

for(l=0;l<NumHidden;l-H-)  { 
net  =  (float)O.O; 
for(j=0  J<Nuinlaputs;j-H-)  { 

net  ■H=  Iiy)utPattem[i][j]*chromosome[k]; 
k-H-; 

} 

net  =  net-chromosome[k]; 
k++; 

hout[l]  =  (float)net/(l+fabs(net)); 

} 

for(l=0;l<NumHidden;l++){ 

out  +=  hout[l]*chromosomePc]; 
k-H-; 

} 

out  =  out  -  chromosome[k]; 
k-H-; 

out  =  (float)out/(Rfiibs(out)); 
printf("output  for  pattern  %d  is  %f\n",i,out); 


} 

/•  CrossOverChromosome,  is  a  disruptive  operator  will  take  two 
induviduals  in  the  population  and  produce  two  children  •/ 


void  CrossOverChromosome(GeneType  •parent  l.GcneType 
•child  l,GeneType  •child2) 

{ 


int  CrossOverPoint;  /•  Randomly  Chosen  Cross  Over  Point  •/ 
int  i;  /•  Counters  •/ 


•parent2,GeneType 
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/*  pick  a  random  point  from  0  to  Number  of  parameter  */ 

srand(seed); 

seed  =  randO; 

CrossOverPoint  =  (int)  (rand()  %  NumberOfParameters); 
/*  Cut  and  Splice  *! 

for(i=0;i<CrossC)verPoint;i-H-)  { 

childl[il  =  parentl[i]; 
cluld2[i]  =  parent2[i]; 

} 

for(i=CrossOverPoint;i<NumberOfParameters;i+-<-){ 
childl[i]  =  parent2[i]; 
child2[ij  =  parentl[ii; 

} 


/•  CopyChromosome  copies  a  given  chromosome  to  another  •/ 

void  CopyChromosome(GeneType  •OldCopy,  GeneType  •NewCopy) 

{ 

inti; 

for(i=0;i<NumberOfParameters;i-H-)  { 

NewCopyp]  =  01dCopy[i]; 

} 


} 


/•  Crossover  will  cross  over  different  chromosomes  selected  with  a 
piobablity  defined  by  the  cross  over  probability  •/ 

void  CrossOver(PopType  *8616016(1,  PopTyp6  *N6w) 

{ 


inti; 

float  probability; 

int  randomChromosome; 

srand(s6ed); 
seed  =  randO; 
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/*  if  the  coin  flip  generates  a  number  less  than  or  equal  to  that  of  the 
cross  over  probability  initiate  crossover,  else  copy  the  induviduals 
into  the  new  population  */ 

for(i=0;i<PopulationSize;i=i+2)  { 

probability  =  ((float)(randO%l  000))/ 1000; 
ii(probability  <=  CrossOverProbability){ 

CrossOverStat-H-; 

CrossC)verChromosome(Selected[i].chromosome,Selected[i+l].chromosome,New[i].chromoso 

me,New[i+l].chromosome); 

) 

else{ 

CopyChr(Mnosome(Selected[i].chromosome,New[i].chromosome); 
CopyChromosome(Selected[i+ 1  ]  .chromosome,New[i+ 1]  .chromosome); 

} 

}/•  end  for  •/ 

/*  if  the  population  size  is  not  even  pass  a  random  chromosome  to  the 
new  population  */ 

if  (PopulationSize%2  !=  0){ 

randomChromosome  =  randO%PopulationSize; 

CopyChromosome(Selected[randomChromosome]  .chromosome,New[randomChromosome]  .c 
hromosome); 

} 

CalculateFitness(New); 


} 


/*  if  the  coin  flip  yeilds  a  probability  less  than  or  equal  to  the 
Mutation  Probability  mutate  gene  */ 

void  MutateChromosome(GeneType  *gene) 

{ 

inti; 

float  probability; 
float  mutateAmount; 
int  power; 

srand(seed); 
seed  =  landQ; 

for(i=0;i<NumberO£Parameters;i++)  { 

probability  =  ((float)(rand0% 1000))/! 000; 
ifl(probability  <  MutateProbability){ 
MutationStat-H-; 


/*  Generates  a  random  mutation  amount  in  the  user  defined  range 
and  either  adds  or  subracts  the  generated  noise  to  the  original 
value  */ 

mutateAmount  =  range*((float)(randO% 1000))/ 1000; 
power  =  rand0%2; 

mutateAmount  =  (float)pow((double)-l,(double)powcr)*mutateAmount; 
gene[i]  =  gene[i]-KGeneType)(inutateAmount); 


} 

}/*endfor*/ 


} 

/*  Mutate  Calls  Mutate  Chromosome  */ 

void  Mutate(PopType  *Population) 

{ 


int  i; 


for(i=0;i<PopulationSi2e;i-H-){ 

MutateChromosome(Population[i].chromosome); 


} 


CalculateFitness(Population); 

} 


/♦  SRSelect  selets  in  one  shot  all  the  orgainsms  for  new  population 
Assumes  that  CalculateFitness  has  already  been  called  */ 

void  SRSelect(PopType  *01d,  PopType  *New) 

{ 


intij; 

float  *sUce;  /*  keeps  the  %  slice  of  the  fitness  for  each  organism  */ 
float  TotalFitness=0;  /*  total  fitness  of  the  populatimi  */ 
float  WheelStop;  /*  Keqjs  the  random  point  where  the  Roulette  Wheel  st(q)s  •/ 
float  delta;  /•  amount  by  which  the  spoke  in  the  vdieel  is  displaced  •/ 
int  picked;  /*  flag  */ 

FTT.F.  /*  Temp  file  variable  */ 


if((^  =  fopen("stat.dat","w+"))  eq  NULL){ 
printfl[''Error:  Cannot  c^ien  file  stat-datW); 
exit(0); 

} 
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/*  Generates  a  random  mutation  amount  in  the  user  defined  range 
and  either  adds  or  subracts  the  generated  noise  to  the  original 
value  */ 

mutateAmount  =  range*((float)(rand0%1000))/1000; 
power  =  rand0%2; 

mutateAmount  =  (float)pow((double)-l,(double)power)*mutateAmount; 
gene[i]  =  gene[i]+(GeneType)(mutateAmount); 


} 

}/*cndfor*/ 

} 

/*  Mutate  Calls  Mutate  Chromosome  V 

void  Mutate(PopType  •Population) 

{ 


inti; 

for(i=0;i<PopulationSize;i-H-)  { 

MutateChromosome(Populationli]. chromosome); 

} 

CalculateFimess(Population), 

} 


/•  SRSelect  selets  in  one  shot  all  the  orgainsms  for  new  population 
Assumes  that  CalculateFitness  has  already  been  called  */ 

void  SRSelect(PopType  ’Old,  PopType  ’New) 

{ 


int  ij; 

float  •slice;  /*  keeps  the  %  slice  of  the  fitness  fo-  each  organism  */ 
float  TotalFitness=0;  /*  total  fitness  of  the  population  */ 
float  WheelStop,  ■'*  Keeps  the  random  point  where  the  Roulette  Wheel  stops  ’/ 
float  delta,  /•  amount  by  which  the  spoke  in  the  wheel  is  displaced  */ 
int  picked;  /•  flag  •/ 

FE^E  •fp;  /•  Temp  file  variable  •/ 


if((fp  =  fopen("stat.dat","w+"))  eq  NIJLL){ 
printfC'Error:  Cannot  open  file  stat.dat\n"), 
exit(O); 

} 
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APPENDIX  C 

OBJECTIVE  FUNCTION  FOR  AUTOMATICALLY  TUNING  THE 

FUZZY  CONTROLLER 


- - *! 

t*  */ 

/♦  The  routines  below  generate  the  control  output  and  compare  */ 
I*  the  root  mean  square  error  of  the  output  produced  by  the  GA  */ 
/*  and  the  actual  ou4)uts.  This  gives  a  measure  of  the  quality 
/*  of  the  fuzzy  membership  function  generated  by  the  G A  */ 

/♦  Sam  Collins,  Dept,  of  Computer  Science,  CAU  */ 

/* - »/ 


#mclude  <math.h> 
/Kinclude  <stdlib.h> 
^include  <string.h> 
^include  cctypc.lP- 
#include  "extemal.h" 


double  evalTriMu(lefl,  center,  right ,  inputValue) 

/•  this  function  returns  the  truth  value  associated  with  the  triangle  •/ 


double 

left; 

double 

center. 

double 

right; 

double 

inputValue; 

{ 

double  slope; 
double  ylnt; 

if((inputValue  >  left)&&(inputValue  center)) 

{ 

slope  =  1 .0/(center  -  left); 
ylnt  =  -slope  •  left; 

}else{ 

slope  =  1.0/(center-right); 
ylnt  =  -slope  •  right; 

} 

retum(slope*inputValue  +  ylnt); 

} 

void  evalInMemMu(member,  inputValue) 

/•  This  function  is  used  to  craluate  the  truth  value  of  a  member.  */ 
struct  membeiObject  •member, 

double  inputValue; 

{ 

switch(member->type) 

{ 

case  TRIANGULAR: 

if(inputValue  =  member->centerVcrtex) 

{ 

member->mu=  1.0; 

}  else  if((inputValue  >  meml^->leftVertex)&& 
(inputValue  <  member->rightVertex)) 

{ 

member->mu  =  evalTriMu(membcr->left  Vertex, 
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} 


member->centerVertex, 

member-->rightVertex, 

inputValue); 

}else{ 

member->mu  =  0.0; 

} 

break; 

case  S_FUNCTION: 

iffinputValue  >=  member->centerVertex) 

{ 

member->mu  =  1.0; 

}else  if{(inputValue  >  member->leftVertex)&& 

(inputValue  <  member->center  Vertex))  { 
member->mu  =  evalTriMu(member->lefl  Vertex, 

member->centerVertex, 

member->centerVertex, 

inputValue); 

}else{ 

member->mu  =  0.0; 

} 

break; 

case  Z_FUNCTION: 

if(input  Value  member^>center  Vertex) 

{ 

member->mu  .0; 

}else  if((input  Value  >  member  ■>centerVertex)&& 

(inputValue  <  member->right  Vertex))  { 
member->mu  =  evalTriMu(member->centerVertex, 

membcr->centerVertex, 

member->rightVertex, 

inputValue); 

}else{ 

member->mu  =  0.0; 

} 

break; 

default: 

§)rintf(stderr,’*\nevallnMemMu:  Error  in  the  case"); 
exit(l); 


void  evalMuO£InVar(C,  Sample) 

f*  this  function  is  used  to  evaluate  input  variables  mu  */ 

struct  controller  ♦C; 

int  Sample; 

{ 

int  i  j; 

i  <  C->numOfInputVar,  1++) 

{ 

foi(j=0;  j  <  C->inputVars[i].numOfMembers;  j++) 

{ 

evalInMemMu(&(C.>inputVars[i].member[j]), 

InputData[Sample][i]); 

} 

} 

} 


void  ZcroMuOK)utVai(C) 

/*  this  function  is  used  to  zero  ou^ut  variables  mu  */ 
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struct  controller  *C; 

{ 

int  ij; 

for(i=0;  i  <  C->numOfDutputVar,  i-H-) 

{ 

for(j=0;j  <C->outputVars[i].niimOfMembers;j*H-) 

{ 

C->outputVars[i].meraber[j].inu=0.0; 

} 

} 

} 


double  mi!iDouble(valuel,  valiie2) 
double  valuel; 
double  value2; 

{ 

if(  valuel  <  valuc2) 

returii(valuel); 

else 


} 


retuni(value2); 


double  maxDouble(  valuel,  value2) 
double  valuel; 
double  value2; 

{ 

if(  valuel  >  value2) 

retum(valuel); 

else 


} 


retuiii(value2); 


double  applyRule(C,R,  mleP,  ruleValue) 

/*  this  fimction  is  used  to  apply  a  rule  if  it  fires  */ 

struct  controller  *C; 

struct  ruleObject  *R; 

struct  ruleComp  ♦•ruleP; 

double  ruleValue; 

{ 

double  truth  Value; 
double  truthForOr, 
truthValue  =  ruleValue; 
while((*ruleP>>type  !=  THEN) 

{ 

switch((*ruleP}->type) 

{ 

case  AND: 


truthValue  =  ininDouble(truthValue, 

C->mputVars[(*ruleP)->varIndex].member[(*ruleP)c>memIndex].mu); 
♦rulcP  =  (♦rul^)^next; 
break; 
case  OR: 


truthForOr  = 

C->inputVars[(  *ruleP)->varIndex]  .member[(  •ruleP)->meinIndex]  .mu; 
♦rulcP  =  (*rul^)->ncxt; 

truthValue  =  maxI)ouble(truthValue,  applyRule(C,R,  ruleP,  truthForOr)); 

/*  recursive  fimction  call  should  point  to  THEN  v/hsa  returns*/ 
if((*ruleP>>type  !=  THEN){ 
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5)rintf(stderT,  "\napplyRuie:  Eiror  in  recursive  call!  ”); 
exit(-l); 

} 

break; 

default: 

fprintf(stderr,’'\napplyRuIe;  Error  in  case”); 
exit(l); 

} 

} 

retuni(truth  Value); 

} 

void  evaIMuOfOutVar(C^) 

/*  This  function  is  utilizes  the  rules  from  the  ruleBase  to  determine  the  */ 

/•  mu  of  the  output  variables  */ 
struct  controller  *0; 
struct  ruleObject  ♦R; 

{ 

int  numOfRules; 

double  ruleValue; 

struct  ruleComp  ‘ruleP, 

for(numOfRules=0;  numOfRules<R->numOfRules;  -HmumOfRules) 

{ 

ruleP  =  &(R->rule[numOfRules]); 
ruleValue  = 

C^inputVars[ruIeP->varlhdex].member[ruleP->memIhdex].mu; 
ruleP  =  ruleP->next; 

ruleValue  =  applyRule(C,  R,  &(ruleP),  ruleValue); 

/•  nileP  should  be  pointing  to  the  T  HEK  component  */ 
if(ruleP->type  !=  THEN) 

{ 

§)rintf(stderr,  "\nevalMuOfDutVar:  Error  should  be  THEN! "); 
exit(-l); 

} 

C->outputVars[ruleP->varIndex].member[ruleP->memIndex]  .mu  = 
maxDouble(rule  Value, 

C->outputVars[ruleP->varIndex]  .member[ruleP->memIndex]  .mu); 

} 

void  calcMemAreaCent(member,  lower,  upper.  Area,  Centroid) 

struct  membeiObject  *member, 

double  lower, 

double  upper, 

double  ♦Area; 

double  *Centroid; 

{ 

double  AreaOne; 
double  AreaTwo; 
double  CentroidOne; 
double  CentroidTwo; 
switch(member->type) 

{ 

case  TRIANGULAR: 

♦Area  =  0.5  *  (member->rightVertex  -  member->leftVertex)  ♦ 
(member->mu); 

CentroidTwo  =((member->rightVertex  -  member->leftVertex)  + 
(member->rightVertex  -  member->centerVertex))/3.0; 
♦Centroid  =  member->right Vertex  -  CentroidTwo; 
break; 

case  S_FUNCTION: 
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AreaOne  =  0.5  *  (member->centerVertex  -  member->ieftVertex)  * 
(member->mu); 

AreaTwo  =  (upper  -  member->centerVertex)  *  (member->mu); 
♦Area  =  AreaOne  +  AreaTwo; 

CentroidOne  =  member->centerVertex  - 

((member->centerVertex  -  member->leftVertex)/3.0); 
CentroidTwo  =  upper  -  (upper  -  member->centerVertex)/2.0; 
if(*Area  >  0.0) 

{ 

♦Centroid  =  ((AreaOne  *  CentroidOne)  + 

(AreaTwo  ♦  CentroidTwo))/  (♦Area); 

}else{ 

♦Centroid  =  0.0; 

} 

break; 

case  Z_FUNCTION: 

AieaOie  =  (  member->centerVertex  -  lower)  ♦  (member->mu); 
AreaTwo  =  0.5  ♦  (meniber->right Vertex  -  member->centerVertex)  ♦ 
(member->mu); 

♦Area  =  AreaOne  +  AreaTwo; 

CentroidTwo  =  nieniber->centerVertex  + 

((member->rightVertex  -  member->centcrVertex)/3.0); 
CentroidOne  =  lower  +  (member*^>centerVertex  -  lower)/2.0; 
if(*Area  >  0.0) 

{ 

♦Centroid  =  ((AreaOne  ♦  CentroidOne)  + 

(AreaTwo  ♦  CentroidTwo))/  (♦Area); 

}else{ 

♦Centroid  =  0.0; 

} 

break; 

defiault: 

fprintf(stderr/\ncalcMeinAreaCent:  Error  in  the  case.\n"); 
exit(l); 

} 

} 

double  CalcOutCentroid(Var) 

/♦  this  function  calculates  the  overall  centroid  of  a  variable  ♦/ 
struct  varObject  ♦Var, 

{ 

int  i; 

double  Area,  Centroid,  TotalArea,  TotalCentArea; 

Area  =  0.0; 

Centroid  =  0.0; 

TotalArea  =  0.0; 

TotalCentArea  =  0.0; 

for(i=0;  i  <  Var->numOfMenibers;  i-H-) 

{ 

calcMeinAreaCent(&(Var->member[i]),  Var->lowerBound, 
Var->upperBound,  &Area,  &Centroid); 

TotalArea  =  TotalArea  +  Area; 

TotalCentArea  =  TotalCentArea  +  (Area  ♦  Centroid); 

i£(TotalArea  >  0.0) 

{ 

retum(  TotalCentArea/TotalArea); 

}else{ 

ietum(0.0); 


} 


} , 


void  getCalcOut(C,  Samp) 

/*  this  function  is  used  to  evaluate  the  calculatef  Out  put  for  a  sample  •/ 
struct  controller  *C; 
int  Samp; 

{ 

inti; 

for(i=0;  i  <  C->numOfDutputVar,  i-H-) 

{ 

CalculatedOutput[Samp][i]  =  CalcOutCentroid(&(C->outputVars[i])); 

} 

} 

double  Zax)CalcOut(numOutVar) 
int  numOutVar, 

{ 

int  numOfSamp^umOfOutputs; 

for(numOfSan:]i>=0;  numOfSamp<NumbeiOfSampIes;  -HuumOfSamp) 

{ 

for(numOfDutputs=0;  numOfOutputs<numOutVar,  -H-numOfOutputs) 

{ 

Ca]culatedOu^ut[niiTrr>f3ai,5p]{jiun3OfDu^uts]=0.0; 

} 

} 

} 

void  calcOutputs(C  Jl) 

/♦  this  function  is  used  to  evaluate  the  fltness  of  the  controller  ♦/ 
struct  controller  *C; 
struct  nileObject  *R; 

{ 

int  San^}; 

ZeroCalcOut(C->numOfOutputVar); 
for(Samp=0;  Samp  <  NumberOfSamples;  Samp+-i-) 

{ 

evalMuOflnVarfC,  Samp); 

ZeroMuOfOutVarfC); 
evalMuOfDutVar(cil); 
getCalcOut(C,  Samp); 

} 

} 

double  ControllerFitness(numOutVar) 

/*  This  function  is  used  to  calculate  the  square  of  the  dififerance  of  */ 

/•  output  data  and  the  calculated  output.  Then  compute  the  Otness.  */ 
int  numOutVar, 

{ 

int  numOfSarap^umOfOutputs; 
intij; 

double  sum=0.0; 
double  mse; 

for(numOfSamp=0;  numOfSan^KNumberOfSamples;  -t-HiumOfSan^)) 

{ 

for(numOfDutputs=0;  numOfOutputs<numOutVar,  ++numOfOutputs) 

{ 

sum  +=((CalculatedOutput[numOfSamp][numOR)utputs]- 
OutputData{numOfSamp][numOfOutputs])* 
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void  getData(nuinInVar,  numOutVar) 

/•  This  function  is  used  to  gel  the  data  from  tlie  Dalai’ ile.  llic  data  will  be 
i*  used  to  calculate  the  ntness  ofUie  indvidual  ♦/ 
uit  nuniInVar, 
int  numOutVar, 

{ 

FILE  *fpGctData; 

int  numOfSamp^umOflnputs^iuxnOfDutputs; 
int  ij; 

ii(numfiles  ==  0) 

fprintffouttp,”  Enter  the  iiajue  of  Uie  datafile. - >  “); 

fscan  f(  mip ,  ”  %s "  ,Da  laF  i  1  e ) , 
if((fpGetData=fopen(DataFile,"r")}  ===  NULL) 

{ 

fprintf(outfp, "Could  not  open  %s\n"4)ataFile); 
cxit(.l); 

} 

fscanftfpGetData,"%d",&NuinberOfSaniples) 

/•  alienate  space  for  two  dimensional  array  InputData  */ 
if{(lnputDala  =  (double  •♦)maJloc 

(Numbcr01Samples*sizeoh;double  ♦)))  =  NULL) 

{ 

fprintf(stderr,"\nNot  enough  memory  for  InputData. \n"); 
cxit(.l); 

} 

for(i  =  0;  i  <  NumberOfSamples;  i-h*) 

{ 

if{(lnpulDala(iJ  =  (double  *) 

malloc(numlnVax*sizeotldouble)))  ==  NULL) 

{ 

fprintf(stdcrT,"\nNot  enough  memory  for  InputData.\n"); 


exit(-l); 


} 


/•  allocate  space  for  two  dimensional  array  OutputData  •/ 
liUOuiputData  =  (double  •*) 

malloc^NumberOfSamples’sizeofCdouble  ^)))  ^  NULL) 

fprintf(stderT,"\nNot  enough  memory  for  OutputData. \n"); 

} 

for(j  =  0;  j  <  NumberOfSamples;  ji-+) 


ill;(OulputJ^atalj)  =  (double  ♦) 


maUocfnumOulVar^sizeotldouble)))  =  NULL) 


fprintf(stdcrr,"\nNot  enough  memory  for  lnputData.\n"); 
exil(-l); 


f*  allocate  space  for  two  dmiensional  array  UulpuU.)aia 
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,  (CalculatedOutput[numOfSamp][nuinOfDutputs]- 

OutputData[numOfSamp][numOfOutputs])); 


mse  =  sqrt(suin/((double)NuinberOfSamples)); 
retum(  1 .0/(  1  .(Hmse)); 


void  getData(nuinInVar,  numOutVar) 

/♦  This  function  is  used  to  get  the  data  fiom  the  DataFile.  The  data  will  be  */ 
/♦  used  to  calculate  the  fitness  of  the  indvidual.  */ 
int  numInVar, 
int  numOutVar, 

{ 

FILE  •§)GetI>ata; 

int  numOfSamp^umOfbiputs^umOfDu^uts; 
intij; 

ih^numfiles  =  0) 

§)rintf(out§),"  Enter  the  name  of  the  datafile. - "); 

fscanfl;inj^,"%s"J)ataFile); 
if((Q)GetData=fopen(DataFile,''r"))  =  NULL) 

{ 

§)rintf(out^,"Could  not  open  %s\n"  JDataFile); 
exit(-l); 

} 

fscanf(fpGetData,  ’*%d",&NumberOfSamples)- 
/*  allocate  space  for  two  dimensional  array  InputData  V 
if((InputData  =  (double  ♦*)malloc 

(NumberOfSamples*si2eof(double  ♦)))  =  NULL) 

{ 

fi)rintf(stderT,"\nNot  enough  memory  for  InputData.Nn"); 
exit(.l); 

} 

for(i  =  0;  i  <  NumberOfSamples;  i-H-) 

{ 

if((lnputData[i]  =  (double  ♦) 


malloc(numInVar*sizeof(double)))  =  NULL) 

^rintf(stderr/\nNot  enough  memory  for  InputData.\n"); 
exit(-l); 


f*  allocate  space  for  two  dimensional  array  OutputData  */ 
if((OutputData  =  (double  •♦) 

malloc(NumberOfSamples*sizeof(double  ♦)))  =  NULL) 

fprintf(stderr,"\nNot  enough  memory  for  OutputE)ata.\n**); 
exit(-l); 

} 

for(j  =  0;  j  <  NumberOfSamples;  j-H-) 

{ 

ifl;(OutputData[j]  =  (double  *) 

malloc(numOutVar*sizeof(double)))  =  NULL) 

§)rmtf(stderr,’*\nNot  enough  memory  for  InputE>ata.\n"); 
exit(-l); 

} 


/*  allocate  space  for  two  dimensional  array  OutputData  ♦/ 
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ifl;(CalculatedOutput  =  (double  **) 

malloc(NumbeiOfSamples*sizeof(double  *)))  =  NULL) 

^)riiitf(stdeir,"\nNot  enough  memory  for  CalculatedOutput.\n"); 
exit(-l); 

for(j  =  0;  j  <  NumberOfSamples;  j++) 

ifr(CalculatedOutput01  =  (double  •)  , 

inalloc(numOutVar*sizeof(double)))  —  NULL) 

^  Q)rmtf(stderr,"\iiNot  enough  memory  for  CaIculatedOutput.\n“); 

exit(-l); 

} 

for(numOfSamp=0;  numO£Samp<NumberOfSamples;  ++numOffiamp) 

^  for(numOfInputs=0;  numOflnputs<JiumInVar,  ++numOfliiputs) 

^  fscajif(fpGetData,"%ir,&InputData[numOfSamp][numOfhq)uts]X 

for(numC)^tputs=0;  numOfOulputs^numOutVar,  ++numOfOu^uts) 

^  fscanh[fpGetData.*%ir,&C)utputData[numOfSamp][numOK>utputs]); 


fclose(fpGetData)y*Close  the  file  •/ 

}/*  End  of  function  getData.  */ 

void  printDatafnumlnVar,  numOutVar) 

/*  This  function  is  used  to  print  the  data  firom  the  DataFile.  */ 
int  numInVar, 
int  numOutVar, 

^  int  numOfSamp^mmOflnputsjiumOfDutputs; 

intij; 

for(numOfSamp=0;  numOfSamp<NumberOfSamples;  -H-numOlSamp) 

{ 

fi)rintf(outQ),"\n"); 

for(numOfInputs=0;  numOflnputs<numinVar,  ++numOfInputs) 

^  ^rintf(outfij."inVar[%2d][%d]  =  %7.31f  ",  numOfSamp+l 

numOflnputs+1,  lnputData[numOfSamp][numOfIiq>uts]); 


) 

^rintf(outfp,"U"); 

for(numOfDutputs=0;  nmnOfDutputs<numOutVar,  -H-numOfOu^uts) 

^  fprintf(outfp,"outVart%2d][%dl  =  %6.31f  ",  numO£SamiH-l, 

mnnOfDu^uts+1,  Ou^tData[numOfSamp][numOfOutputs]); 

} 

} 

} 

int  gethiVarIndex(C,  name,  num) 

Struct  ccmtroller  *C\ 
char  •name; 
int  num; 

{ 

inti; 
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for(i=0;  i  <  C->numOfInputVar,  i-H-) 

{ 

if(!strcmp(name,  C->inputVars[i].varName)) 
retum(i); 

} 

5)rintf(stderT,  "\ngetInVarIndex:  Error  in  rule[%d]l  "^um); 

^rintf(stderr,  "Variable  *%s'  does  not  exist!\n\n",  name); 
exit(-l); 

} 

int  getInMemIndex(C,  name,  num,  varindex) 

struct  controller  ’C; 

char  *name; 

int  num; 

int  varLadex; 

{ 

inti; 

foi(i=0;  i  <  C->ii^)utVars[varIndex].numOfMembers;  i-H-) 

{ 

if(!strcmp(name,  C->ii^)utVars[varIndex],member[i].memName)) 
retum(i); 

} 

Q)rintf(stderr,  ’’\ngetInMemIndex:  Error  in  nile[%d]!  "^lum); 

Q}rintf(stdeiT,  "Member  *%s'  does  not  exist!\n\n",  name); 
exit(-l); 

} 

int  getOutVarIndex(C,  name,  num) 
struct  controller  *C; 
char  ♦name; 
int  num; 

{ 

inti; 

for(i=0;  i  <  C->numOfOutputVar,  i-H-) 

{ 

if(lstrcmp(name,  C->outputVars[i].varName)) 
retum(i); 

} 

fprintf(stderT,  "\ngetOutVarIndex:  Error  in  rule[%d)!  "^lum); 
fi>rintf(stderr,  "Variable  ^s'  does  not  exist!\n\n",  name); 
exit(-l); 

} 

int  getOutMemIndex(C,  name,  num,  varindex) 

struct  controller  ♦€; 

char  ♦name; 

int  num; 

int  varindex; 

{ 

inti; 

for(i^;  i  <  C->outputVars[varIndex].numOfMembers;  i-H-) 

{ 

if(!strcmp(iiame,  C->outputVai8[varIndex].memberfi].memName)) 
reti]m(i); 

} 

Q)rintfl[stdeir,  ”\nget0utMeinlndex:  Error  in  iule[%d]!  "jium); 

Q>rintf(stderr,  'Member  '%s'  does  not  exist!\n\n",  name); 
exit(.l); 

} 
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int  buildRuleBase(C4iileP,fpGR^mn,symbol) 
struct  controUer  •€; 
struct  ruleComp  ♦ruleP; 

¥RE**fpGR, 
intnum; 
char  *symbol; 

{ 

char  symbolIS[MAX]; 
if{!strcmp(symbol,  "ff*")) 

{ 

ruieP->type  =  ff ; 

fscanf(*^GR,’’%s",  ruleP->varName); 

ruleP->varIndex  =  getInVarIndex(C,  ruleP->varName,  num); 

fscanfi:*^R,"%s",  symboIIS); 

if(strcmp(symbolIS,  "IS")) 

{ 

§)rintf(stderr,  "\nbuildRuIeBase:  Error  in  rule  %d"^um); 
^rmtf(stderr, ".  Expecting  IS!\n"); 
cxit(.l); 

} 

fscanf(*§)GR,"%s",  ruleP->memNanie); 
ruleP*>n^mIndex=getIriMemIndex(C,  iuleP->meinName , 

num,  ruleP->varIndex); 


retum(l); 

}else  if(!strcmp(symbol,  "AND")){ 
ruleP->^pe  =  AND; 
fscanf(*fi^R,"%s",  ruleP->varName); 
ruleP->varIadex  =  getInV8rIndex(C,  ruleP->varName,  num); 
fscanft*fpGR,"%s",  symboIIS); 
if(strcmp(symbolIS,  "IS")) 

{ 

fprintf(stderr,  "\nbuildRuleBase:  Error  in  rule  %d"4nim); 
fj^tf(stderr, ".  Expecting  IS!\n"); 
exit(-l); 

} 

fscanf{*§)GR,"%s",  ruleP“>memName); 
nileP->meniIndex=getIiiMemIndex(C,  ruleP->memName , 

num,  ruleP->varIndex); 


retum(l); 

}else  if{!strcII^)(symbol,  "OR")){ 
ruleP->type  =  OR; 

fscanfr*fi^R,"%s",  ruleP->varName); 

ruleP->varIndex  =  getInVarIndex(C,  ruleP->varName,  num); 

fscanf(*fpGR,"%s",  symboIIS); 

if(strcmp(symbolIS,  "IS")) 

{ 

Q)rintf(stderr,  "VnbuildRuleBase:  Error  in  rule  %d"4ium); 
fl^tf(stderr, ".  Expecting  IS!\n"); 
cxit(-l); 

} 

fscan5*fpGR,"%s",  ruleP->memName); 
ruleP->memIndex=getInMemIndex(C,  ruleP->mcmName , 

num,  rulcP«>varIndex); 


retuin(l); 

}else  if(!strcn^)(symbol,  "THEN")){ 
rulcP->type  =  THEN; 
fscanf(*fi^R,"%s",  ruleP->varName); 
ruleP->varIndex  =  getOutVarIndex(C,  ruleP->varName,  num); 
fscanf(*ft)GR,"%s",  symboIIS); 
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if(strcmp(symbolIS,  "IS")) 

{ 

fprintf(stderr,  "\nbuildRuleBase:  Error  in  rule  %d"4ium); 
f^rmtf(stderr, Expecting  IS!\n"); 
exit(-l); 

} 

fscanf(*§)GR,"%s",  ruleP->meinName); 
ruleP->memIndex  =  getOutMemIndex(C,  nileP->meniName, 

num,  ruleP->varIndex); 

retum(0); 

}else{ 

fprintf(stderr,  "\nbuildRuIeBase:  Error  in  rule  %d  file. ",  num); 
fprintf(stderr,  "Check  fommt!"); 
exit(-l); 

} 

} 

struct  ruleComp  *allocateCon^)() 

{ 

struct  ruleCon^  ‘RC; 
if((RC  =  (struct  ruleCon^)  *) 

iiialloc(sizeof(struct  ruleComp)))=NULL) 

{ 

fprintf(stderT,"\nallocateComp:  not  enough  memory  for  mles.\n"); 
exit(-l); 

} 

rctum(RC); 

} 

getRuies(C,R) 

/•  This  function  is  used  to  get  the  rule  fium  RuleFile.  */ 
struct  controller 
struct  ruIeObject  ♦R; 

{ 


FILE 

•fpGetRules; 

int 

numOfRules; 

char 

symbol[MAX]; 

int 

endFlag; 

int 

ne\\^ode; 

struct  ruleComp  *P; 
if(numfiles  =  0)  { 

fprintf(outfp,"  Enter  the  name  of  the  rule  file. "); 
fprintf(outfp," - > "); 

} 

fscanf(infp,"%s"3^uleFile); 
ifl;(fpGetRules=fopen(RuleFile,"r"))  =  NULL) 

{ 

§)rintf(stdeiT,"Could  not  open  %s\n",RuleFile); 
cxit(-l); 

} 

fscanf(^)GetRules,"%d",&(R->numOfRules)); 
if((R->iule  =  (struct  ruleComp  •) 

malloc((R->numOfRules)*si2eof(struct  ruleComp)))=NULL) 

{ 

fprintf(stderr,"\ngetRules:  not  enough  memory  for  niles.Nn"); 
exit(-l); 

} 

for(numOfRules=0;  numOfRules<R->numC)fRules;  -KnumOfRules) 

{ 

endFlag  =  0; 
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newNode  =  0; 

P  =  &(R->rule[numOfRules]); 

!feof(fpGetRules)  &&  lendFlag) 

{ 

fscanfl^QjGetRules/yos”,  symbol); 
if(!strcmp(symbol,  ’'END"^ 

{ 

endFlag  =  1; 

}else{ 

newNode=buildRuleBase(CJP,&fpGetRuIes, 
numOfRules,  symbol); 

if(newNode) 

{ 

P->next=allocateComp(); 

P=P->next; 

}else{ 

P->next=NULL; 

} 

} 

} 

prmtf("\n"); 

} 

fclose(fpGetRuIes)/*Close  the  file  •/ 

}/♦  End  of  function  getRules.  ♦/ 


void  printRules(R) 
struct  ruleObject  *R; 

{ 

struct  ruleComp  *P; 

int  numOfRules; 

4)rintf(stdout,’’\n\n’'); 

foifnumOfRules=0;  numOfRules<R->numOfRules;  ++numOfRules) 

{ 

P  =  &(R->rule[numOfRules]); 
fprintf(stdout/RULE[%d]  =  "^lumOfRules); 
vAule  (P!=NULL) 

{ 

switcli(P->type) 

{ 

case  IF: 

fprintf(stdout,"IF  "); 
break; 
case  AND: 

Q)rintf(stdout,"AND "); 
break; 

case  OR: 

§)rintf(stdout,"OR "); 
break; 
case  THEN: 

Q)rintf(stdout,"THEN "); 
break; 

default 

fprintf(stderr,"\aprintRules: "); 

fprintf(stderr,’*In  case.\n"); 

exit(-l); 

} 

fprintf(stdout,"%s  IS  %s ",  P->varName,  P->memNamc); 
P=P->next; 

} 
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Q)rintf(stdout,  '’END\n" ); 


} 

§)rintf(stdout,”\n\n’’); 

} 

void  freeRuIes(R) 
struct  ruleObject  ^R; 

{ 

struct  nileComp  *P; 
struct  ruleComp  ♦tmp; 
int  numOfRules; 

for(numOfRules=0;  numOfRules<R->numOfRules;  ++numO£Rules) 

{ 

P  =  &(R->niIe[numOfRules]); 
while  (P->next!=NULL) 

{ 

tmp  =  P->next; 

P->next=P->next->next; 

fiee(tmp); 

} 

} 

free(R->ruIe); 

} 


void  evalMembers(memberOne,  memberTwo) 

/•  This  function  is  used  to  evaluate  the  validity  of  two  member  functions  •/ 
struct  memberObject  •memberOne; 
struct  memberObject  *memberTwo; 

{ 

switcb(membeiOne->type) 

{ 

case  TRIANGULAR: 

switch(memberTwo->type) 

{ 

case  TRIANGULAR: 
break; 

case  S_FUNCTION: 

fprintf(stderr,"\nevalMembers:  Error  invalid "); 
fprintf(stderr»"combination  of  member  functions!"); 
fi)rintf(stderr,"\n\n\tS_FUNCTION ...  TRIANGULAR  \n"); 
exit(l); 
break; 

case  Z_FUNCTION: 
break; 

default: 

§)rintf(stderr,"\nevalMembers:  Error  in  inner  case"); 
exit(lX 

} 

break; 

case  S^FUNCnON: 

switch(memberTwo->type) 

{ 

case  TRIANGULAR: 
break; 

case  S^FUNCTION: 

fprintf(stderr,"\nevalMembers:  Error  invalid  "); 
fprmtf(stderr,"combination  of  member  functions!"); 
^tf(stderr,"\n\n\tS^FUNCT[ON ...  S_FUNCT[ON\n"); 


} 

} 


exit(l); 

breaJc; 

case  Z_FUNCTION: 
break; 

default: 

Q)rintf{stderr,"\nevalMembers:  Error  in  inner  case”); 
exit(l); 

} 

break; 

case  Z_FUNCTION: 

switch(memberTwo->type) 

{ 

case  TRIANGULAR: 

§)rmtf(stdeiT,"\nevaIMembers:  Error  invalid "); 
Qjrintffstderr/combination  of  member  fimctionsr); 
fprintf(stderr/\n\n\tTRIANGULAR ...  Z^FUNCTlONVn”); 
exit(l); 
break; 

case  S_FUNCTION: 

fprintf(stderr,"\nevalMembers:  Error  invalid  ”); 
Q)rintf(stderr,"combmation  of  member  fimctionsr); 
§)rintf(stden*,"\n\n\tS_.FUNCTION ...  Z  FUNCTION^"); 
exit(l); 
break; 

caseZ_FUNCTION: 

§)rintf(stderr,"\nevalMembers.  Error  invalid  ”); 
§)rintf(stderr, "combination  of  member  fimctions!"); 
fprmtf(stderr,"\n\n\tZ^FUNCT[ON ...  Z^FUNCTION\n"); 
exit(l); 
break; 

de&ult: 

fprintf(stderr,"\nevalMembers:  Error  in  iimer  case"); 
cxit(l); 

} 

break; 

de&ult 

fprintf(stderr,"\nevalMembers:  Error  in  outer  case"); 
exit(l); 


void  checkMem(TheVar) 

/•  This  fimction  will  be  used  to  determine  the  number  the  of  criteria  used  */ 
/*  to  calculate  the  fitness  V 
struct  varObject  *TheVar, 

{ 

intjjc; 

for(j=0;  j  <  ThcVar->numOfMembers;  j++) 

{ 

for(k=^-l;  k>=0;  k-) 

{ 

evalMembers(TheVar->membcrlj], 

TheVar->member[k]); 

} 

} 

} 

int  getType(nameM) 
diar  •nameM; 

{ 
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int  choice; 


/*  varible  that  holds  the  *f 

f*  choice  of  the  user. 


/* 

*  This  part  of  the  program  will  print  the  list  of  choices  for  the  user. 

V 

if(numfiles  =  0){ 

/*^rint^outfp,"\n");V 

fprintf(outfp,"\tMember  function  *%s'  is  of  what  type:\n\n", 
nameM); 

fprintf(outfp,"\t\tl.  Triangular.\n’’); 

^rintftout4),"\t\t2.  S-Function,\n"); 

§)rintfl;out^/\t\t3.  Z-Function.\n’’); 

Qxrintf(outfi),"\n'’); 

§>rintf(outfi),’’\tEnter  your  choice: "); 

} 

fscanf(infp/%d",&choice); 
if(nuinflles  ==  0) 

^rint£(outfp,"\n"); 

/* 

*  This  part  of  the  program  will  return  the  value  for  the  function 

*  that  corresponds  to  die  users  choice. 

V 

if  (choice  =  1 ){ 

ietuin(TRIANGULAR);  /♦  Triangular  •/ 

}  else  if  (choice  =  2){ 

retuni(S_FUNCTION);  /*  S-Function  */ 

}else  if  (choice  —  3){ 

retum(Z_FTJNCTION);  /*  Z-Function  */ 

}else  { 

fprintf(stdeiT,"\nIncorrect  entry  for  type!"); 
exit(.l); 

}  } 

void  allocateInputVar(C) 

/*  This  function  is  used  to  allocate  space  for  the  input  variables  */ 
struct  controller  *C\ 

{ 

int  ij;  /*  for  loop  variables  */ 

/*  get  the  number  of  input  variables  ♦/ 
if(numfiles  =  0){ 

^rintf(outfp,"\n  Enter  number  of  input  variables  "); 
fprintf(outfp," - > "); 

} 

fscanf(infp,"%d",&(C->numOfInputVar)); 
if((C->inputVars  =  (struct  varObject  ♦) 

malloc((C->numOfInputVar)*sizeof(struct  varObject)))=NULL) 

{ 

§>rintf(stderr,"not  enough  memory  for  %s.\n","inputVars"); 
exit(-l); 

] 

for(i=0;  i  <  C->numOfInputVar, !++) 
if(numfiles  =  0){ 

fprmtf(outfp,"\n  Enter  name  of  input  variable"); 
fprintf(out^,"[%d] ",  i+1); 
fi)rintf(outfp," - > "); 

} 

fscanf(infp,"%s",C->inputVars[i].varName); 
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} 


if^numfiles  =  0){ 

fprintf(outfp,’’\n  Enter  the  lower  bound  ”); 
fprintf(outfp/for  the  variable "); 
fj5rintf(outfj),’"%s'  ->  ”,  C->inputVars[i].varName); 

} 

fscanf(inft>,”%ir,&(C->inputVars[i] .  lowerBound)); 
if(nun^es  =  0){ 

fj^tf(out§),"\n  Enter  the  upper  bound  ”); 
fprintf(outfp,'’for  the  variable  ”); 
fprintf(outfp/'%s'  ->  ”,  C->inputVars[i].varName); 

} 

fscanf(infp,”%ir,&(C->inputVars[i].upperBound)); 

if(  C->mputVars[i].upperBound  <  C->inputVars[i]. lowerBound) 

{ 

fprintf(stderr,"\n  Upper  bound  is  greater  than "); 
fprintf(stderr, "lower  bound  of  variable  "); 
^rintf(stderr,"'%s',\n",  C->inputVars[i].varName); 
exit(-l); 

} 

C->inputVars[i].delta  =  C->inputVars[i].upperBound  - 
C->inputVars[i],lowerBound; 

ifl^numfiles  ==  0){ 

fprintf(out^,"\n  Enter  the  number  of  bits "); 
fprintf(outfi),"for  the  variable "); 

$)rintfl;outfi),"'%s'  -> ",  C->inputVars[i].varName); 

} 

fscanf(infp,"%d",&(C->inputVars[i].numBits)); 
if(numfiles  =  0){ 

fprintf(outip,"\n  Enter  number  of  member  functions "); 

^rintf(out^,"for  variable "); 

fprintf(outfi),"'%s'  ~> ",  C->inputVars[i].varName); 

} 

fscanf(infp,"%d",&(C->inputVars[i]  .numOfMembers)); 

/♦  allocate  space  for  member  functions  for  the  variable  •/ 
if((C->inputVars[i].member  =  (struct  memberObject  •) 
malloc{(C->inputVars[i]. numOfMembers)  * 
sizeof(struct  membeiObject)))=NULL) 

{ 

fprintf(stderr,"not  enough  memory  for  the  "); 
fprintf(stderr,"member  functions  of  variable  "); 
fprintf(slderr,"’%s'.\n",  C->iiiputVars[i].varName); 
exit(-l); 

} 

if(numfiles  =  0) 
fprintf(outfp,"\n"); 

/♦  get  name  and  type  for  member  functions  */ 
for(j=0;  j  <  C->inputVars[i]. numOfMembers;  j-H-) 

{ 

if(nurafiles  =  0){ 

fprintf(outQ),"\tEnter  name  of  member "); 

fprmtf(outfp, "function  %2d  ",  j+1); 
fprintf(outfj>,"  -> "); 

} 

fscanf(infp,"%s"  ,C->inputVars[i]  .memberlj]  .memName); 
C->inputVars[i].member[j].type  = 

getType(C->inputVars[i].member|j].mcmName); 

} 

/•  check  for  invalid  combination  of  member  flmctions  •/ 
checkMem(C->inputVars[i]); 
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} 


void  freeInputVar(C) 

/*  This  function  is  used  to  deallocate  space  for  the  input  variables  */ 
struct  controller  *C; 

{ 

inti; 

for(i=0;  i  <  C->numOflnputVar,  i-H-) 

{ 

free(C->inputVars[i].member); 

} 

free(C->inputVars); 

} 

void  allocateOutputVar(C) 

/*  This  function  is  used  to  allocate  space  for  the  ou^ut  variables  */ 
struct  controller  *C; 

{ 

int  i  j;  /*  for  loop  variables  */ 

/♦  get  the  number  of  output  variables  */ 
if(numfiles  =  0){ 

fprintf(outQ)/\n  Enter  number  of  output  variables  "); 

^rintf(outfp," - >  "); 

} 

fscanf(infp,"%d’’,&(C->numOfOutputVar^; , 
if((C->outputVars  =  (struct  varObject  *) 

malloc((C->numOfDutputVar)*sizeof(struct  varObject))>=NULL) 

{ 

fprintf(stdeiT,"not  enough  memory  for  %s.\n" /output Vars"); 
exit(-l); 

} 

for(i=0;  i  <  C->numOfOutputVar,  i++) 

{ 

if(numfiles  =  0){ 

^rintf(outfp/\n  Enter  name  of  output  variable”); 
fprintf(outfp/[%d] ",  i+l); 
f^rintf(outfp," - >  "); 

} 

fscanf(infp,"%s",C->outputVars[i].varName); 
if(numfiles  =  0){ 

fprintfroutfp,"\n  Enter  the  lower  bound  "); 
fprintf(outf^/for  the  variable "); 
fprintf(outf^,"'%s'  — >  ",  C->outoutVars[i].varName); 

} 

fscanf(  infp,  "%lf’  ,&(C->output  Vars[i] .  lowerBound)); 
if(numfiles  =  0){ 

fprintf(outfp,"\n  Enter  the  upper  bound  "); 
fprintf(outfj),"for  the  variable "); 
f^rintf(outf^,"*%s'  — >  ",  C->outputVars[i].varName); 

} 

fscanf(in^,"%ir,&(C->outputVars[i]  .upperBound)); 

if(  C->outputVars[i].upper^und  <  C->outputVars[i].lowerBound) 

{ 

fprintf(stderr,"\n  Upper  bound  is  greater  than "); 
fi)rintf(stderr,"lower  bound  of  variable  "); 
fprintf(stdeiT,"‘%s'.\n",  Cc>outputVars[i].varName); 
cxit(-l); 

} 

C«>outputVars[i].delta  =  C->outputVars[i].upperBoxmd  - 
C->outputVars[i].  lowerBound; 
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if(nuiiifiles  ==  0){ 

fprintf(outfp,”\n  Enter  the  number  of  bits  ”); 
fi)rintf(outfp/for  the  variable  ”); 
fprintf(outfj),"'%s'  -> ",  C->outputVars[i].varName); 

} 

fscanf(inj5),'*%d",&(C->outputVars[i].numBits)); 
if(numfiles  =  0){ 

^rintf(outfp,"\n  Enter  number  of  member  functions  "); 

^rintf(outfp,"for  variable "); 

fi)rintf(outfp,"'%s'  -> ",  C->outputVars[i].varName); 

} 

fscanf(in§),"%d",&(C->outputVars[i].niimOfMembers)); 

/*  allocate  space  for  member  functions  for  the  variable  */ 
if((C->outputVars[i].member  =  (struct  memberObject  •) 
malloc((C->outputVars[i].numOfMemb^)  * 
si2eof(struct  memberObject)))=NULL) 

{ 

Q)rintf(stderr,"not  enough  memory  for  the  "); 
fprintf(stderT, "member  functions  of  variable  "); 
§)rintf(stderr,"'%s’.\n",C->outputVars[i].  varName); 
exit(-l); 


if(numfiies  =  0) 

Q>rintf(outQ),"\n"); 

/♦  get  name  and  type  for  member  functions  •/ 
for(j=0;  j  <  C->outputVars[i].numOfMembers;  j++) 

{ 

if(numfiles  =  0){ 

fprintf(outfp,"\tEnter  name  of  member "); 

fprintf(outfp,"function  %2d ",  j+l); 
fprintf(out^,"  -> "); 

} 

fscanf(infp,"%s",  C->outputVars[i].memberlj].memName); 
C->ou^utVars[i].memberU].type  = 

getTyp^C->ou^utVars[i].member[j].memName); 

} 

/*  check  for  invalid  combination  of  member  fimctions  */ 
checkMem(C->outputVars[i]); 

} 

} 

void  freeOutputVar(C) 

/*  This  function  is  used  to  deallocate  space  for  the  output  variables  */ 
struct  controller  *C; 

{ 

inti; 

for(i=0;  i  <  C->numOfDutputVar,  i++) 

{ 

free(C->outputVars[i].member); 

} 

free(C->outputVars); 

} 

void  defineController(C) 

/*  This  function  is  used  to  allocate  space  for  the  variables  */ 
struct  controller  *C; 

{ 

i^numfiles  =  0){ 

fprmtf(outfp,"\n  Enter  name  of  controller "); 

^rintf(outfp," - "); 


} 

fscanf(infp,"%s",C->controllerName); 

allocateInputVar(C); 

allocateOutputVar(C); 

} 


void  freeVariables(C) 

/•  This  function  is  used  to  allocate  space  for  the  variables  */ 
struct  controller  *C; 

{ 

freelnputVarCC); 

fieeOutputVar(C); 

} 


void  printMemType(Type) 
intType; 

{ 

switch(Type) 

{ 

case  TRIANGULAR: 

§nintf(outQ),"\n\tTiangular "); 
break; 

case  S_FUNCTION: 

5>rintf(outQ),'’\n\tS-Function:"); 

break; 

case  Z_FUNCTION: 

5)rintf(outQ),  "\n\tZ-Function:" ); 
break; 

default: 

Q)rintf(stderr,"\nprintMemType:  Error  in  the  case"); 
exit(l); 

} 

} 

void  printValues(member) 

/•  This  function  is  used  to  print  the  values  of  the  verticies  */ 
struct  memberObject  •member, 

{ 

fprintf(outfp,"\n"); 

switch(member->type) 

{ 


case  TRIANGULAR: 

Q)rintf(outQ),"\n\t\t\t\tleftVertex  =%.3ir, 

member->leftVertex); 

§)rmtf(outfp,"\n\t\t\t\tcenterVertex  =  %.3ir, 

member->centerVertex); 

fpiintf(outfp,"\n\t\t\t\trightVei1^  =%.3ir, 

member->rightVertex); 

break; 

caseS^FUNCTION: 

5)rintf(out^,"\n\t\t\t\tleftVertex  =%.3ir, 

member->leftVertex); 

§MTntf(out§),"\n\t\t\t\tcenterVertex  =  %.3ir, 

member-ocenterVertex); 

break; 

case  Z^FUNCTTON: 

§)rintf(outfp,"\n\t\t\t\tcenterVertex  =  %.3ir. 
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} 


default: 


} 

fprintf(outfp,"\n"); 


member->centerVertex); 
^rmtf(outfp,"\n\t\t\t\trightVertex  =  %,3ir, 

inember->rightVertex); 


break; 


5)rintf(stderr,"\nprintValues:  Error  in  the  case"); 
exit(l); 


void  printInVars(C,  pValue) 

/*  this  function  is  usai  to  check  the  input  variables  ♦/ 
struct  controller  *C; 
int  pValue; 

{ 

int  i  j; 

QjrintiR;outfp,"\n\n - INPUTS - \n"); 

§)rintf(outfp,"\nnumber  of  input  varibles  =  %d",  C->numOfInputVar); 
foi(i=0;  i  <  C->numOfIi^)utVar,  i++) 

{ 

5)rintf(outQ),"\n\nvariable[%d]  =  '%s'  and  it  has  ", 
i+l,  C->inputVars[i].varName); 

§)rintf(out^,"%d  members:",  C->inputVars[i].numOfMembers); 
fprintf(out^,"\n"); 

for(j=0;  j  <  C->inputVars[i].numOfMembers;  j-H-) 

{ 

printMemType(C->inputVars[i]  .memberjjJ.type); 
fprintf(outfp,"\tmemb^%d]  = ",  j+1); 
fprintf(out§),"%s "  ,C->inputVars[i]  .member[j  ] .  memName); 
if(pValue=l) 

printValues(C->inputVars[i],member|j  ]); 

} 

fprintf(outfp,"\n\nlower  bounds  =  %ir, 

C->inputVars[i].  lowerBound); 

fprintf(outfp,"\nupper  bounds  =  %ir, 

C->ii^)utVars[i].upperBound); 

fprintf(outfp,"\ndelta  =  %!!", 

C->inputVars[i].delta); 

fjprintf(outfp,"\nnum.  of  bits  =  %d", 

C->inputVars[i].numBits); 


} 

§)rintf(outfp,"\n\n"); 


void  printOutVars(C,  pValue) 

/•  this  fimction  is  used  to  check  the  ou^ut  variables  •/ 
struct  controller  *C; 
int  pValue; 

{ 

int  i  j; 

§)rmtf(outQ),"\n\n - OUTPUTS - \n"); 

fprintf(outQ),"\nnumber  of  output  varibles  =  %d",  C->nximOfOutputVar); 
for(i=0;  i  <  C->numOfDutputVar,  1++) 

{ 

§)rintf(outfp,"\n\nvariable[%d]  =  and  it  has ", 

i+1,  C->outputVars[i].varName); 
fprintf(outfp,"%d  members:",  C->outputVars[i].numOfMembers); 
fprintf(outQ),"\n"); 

for(j^;  j  <  C->outputVars[i].numOfMembers;  j++) 
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{ 

printMemType(C->outputVars[i].member[j].type); 
fprintf(oulip,"\tmember[%d]  =  ”,  j+1); 

^rmtf(out^,'’%s  ",C->outputVars[i].member|j].meniName); 
if(pValue=l) 

printValues(C->outputVars[i].member|j]); 

} 

fprintf(outfp,’^n\nlower  bounds  =  %ir, 

C->outputVars[i].lowerBound); 

fprint£(outfp,"\nupper  bounds  =  %ir, 

C->outputVars[i].upperBound); 

5)rintf(outfp,"\ndelta  =  %ir, 

C->outputVars[i].delta); 

§>rintf(outfp/\nnum.  of  bits  =  %d", 

C->outputVars[i].numBits); 

} 

§>rintf(outfp,"\n\n"); 


void  print Variables(C  ,  P) 

/•  this  function  is  used  to  print  all  the  variables  */ 

struct  controller 

intP, 

{ 

Q5rintf(outfp,"\n\nUCONTROLLER: "); 

^rintf^outf^/Vos",  C->controllerName); 
printInVars(C,  P); 
printOutVars(C,  P); 

} 

void  evalMembeiCrite(member,  count,  valid) 

/•  This  function  is  used  to  evaluate  the  criteria  for  a  single  member  */ 
struct  membeiObject  *member, 
double  *count; 
double  ’valid; 

{ 

switch(member->type) 

{ 

case  TRIANGULAR: 

’count  =  ’count  +  3.0; 

ifi;member->centerVertex  >=  member->iefl  Vertex) 
•valid  =  ’valid  +  1 .0; 

if(member->right  Vertex  >  member->left  Vertex) 

’valid  =  ’valid  +  1 .0; 

if(member->right  Vertex  >=  member->centerVertex) 
•valid  = ’valid  +  1.0; 

break; 

case  S^FUNCTtON: 

•count  =  ’count  +  1 .0; 

if(member->center  Vertex  >=  member->leftVertex) 
•valid  =  ’valid  +  1.0; 

break; 

caseZ_FUNCT10N: 

•count  =  ’count  +  1,0; 

if(member->rightVertex  >=  member->cenlerVertex) 
♦valid  =  ’valid  +  1.0; 

break; 

default: 

fprintf(stderr,"\nevalMembeiCrite:  Error  in  the  case"); 


exit(l); 


void  coinpareMemberCrite(memberOne,  memberTwo,  count,  valid) 

/•  This  function  is  used  to  evaluate  the  criteria  of  two  member  */ 

struct  memberObject  •memberOne; 

struct  memberObject  *memberTwo; 

double  *count; 

double  *valid; 

{ 

switch(memberOne->type) 

{ 

case  TRIANGULAR: 

switch(memberTwo^type) 

{ 

case  TRIANGULAR: 

*coimt  =  •count  +  6.0; 

if(memberOne-’>left  Vertex  >=  memberTwD->leflVertex) 
•valid  =  •valid  +  1 .0; 

if(memberOne->centerVertex  >  memberTwo->leftVertex) 
•valid  =  •valid  +  1 .0; 

if(membeiOne->centerVertex  >  memberTwo*>centerVertex) 
•valid  =  •valid  +  1 .0; 

if(memberOne->rightVertex  >  memberTwo->leflVertex) 
•valid  =  •valid  +  1 ,0; 

if(memberOne->right  Vertex  >  memberTwo->centerVertex) 
•valid  =  •valid  +  1.0; 

if(memberOne->right  Vertex  >=  memberTwo->rightVertex) 
•valid  =  •valid  +  1.0; 


break; 

case  S_FUNCTION: 

fprintf{stderr/\ncompareMemberCrite:  Error  invalid "); 
fjprintf(stderr,"combination  of  member  functions!**); 
fprintf(stderr,’’\n\n\tS_FUNCT[ON ...  TRIANGULAR  \n-); 
exit(l); 
break; 

case  Z_FUNCTION: 

•count  =  •count  +  3.0; 

if(meraberOne->centerVertex  >  memberTwo->centerVertex) 
•valid  =  *valid  +  1.0; 

if(memberOne->rightVertex  >  memberTwo->centerVertex) 
•valid  =  *  valid  +  l.O; 

if(memberOne->right Vertex  >=  memberTwo->rightVertex) 
•valid  =  •valid  +  1.0; 

break; 

default: 

§)rintf(stderT,'*\ncompareMembeiCrite:  Error  in  inner  case"); 
exit(l); 

} 

break; 

case  S_FUNCTION: 

switch(memberTwo->type) 

{ 

case  TRIANGULAR: 

•count  =  •count  +  3,0; 

if(membeiOne->leflVertex  >=  memberTw<»leftVertex) 
•valid  =  •valid  +  1.0; 

if{memberOne->centerVertex  >  memberTwo->leftVertex) 
•valid  =  •valid  +  1.0; 
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} 

} 


if(memberOne->centerVertex  >  memberTwo->centerVertex) 
♦valid  =  ♦valid  +  1.0; 

break; 

case  S_FUNCTION: 

fprintj[{stderr,"\ncompareMemberCrite:  Error  invalid "); 
fprintf{stderr,"combination  of  member  functions!"); 
$rintf(stderr,"\n\n\tS^FUNCTION ...  S^FUNCTION\n"); 
exit(l); 
break; 

case  Z_FUNCTION: 

♦count  =  ♦count  +  1 .0; 

if(membeiOne->centerVertex  >  memberTwo->centerVertex) 
♦valid  =  ♦valid  +  1.0; 

break; 

detimlt: 

Q)rintf(stderr,"\ncompareMemberCrite:  Error  in  inner  case"); 
exit(l); 

} 

break; 

caseZ_FUNCTION: 

switch(memberTwo->type) 

{ 

case  TRIANGULAR: 

§)rintf{stderr,"\ncompareMemberCrite:  Error  invalid  ”); 
fi)rintf(stderr,"combination  of  member  functions!"); 
^iintfl:stdeiT,"\n\n\tTRIANGULAR ...  Z_FUNCTION\n"); 
exit(l); 
break; 

case  S_FUNCTION: 

fprintf(stderr,"\ncompareMemberCrite:  Error  invalid "); 
fprmtf(stderr,  "combination  of  member  functions!"); 
5)rmtf(stderr,"\n\n\tS__FUNCTO  ...  Z^FUNC'nON\n"); 
exit(l); 
break; 

caseZ^FUNCnON: 

fprintf(stderr,"\ncompareMemberCrite:  Error  invalid "); 
fi)rintf(stderr,"combination  of  member  functions!"); 
fi}rintf(stderr,"\n\n\tZ_FUNCT[ON ...  Z_^FUNCTION\n"); 
exitd); 
break; 

default: 

fprintf(stderr,"\ncompareMemberCrite:  Error  in  inner  case"); 
exit(l); 

} 

break; 

default: 

Q)rintf(stderr,"\ncompareMemberCrite:  Error  in  outer  case"); 
exit(l); 


evaluateCrite(TheVar,  count,  valid) 

/♦  This  function  will  be  used  to  determine  the  number  the  of  criteria  used  ♦/ 

/♦  to  calculate  the  fitness  */ 

stnict  varObject  ♦TheVar, 

double  ♦count; 

double  ♦valid; 

{ 

intjjc; 

for(j=0;  j  <  TheVar->numOfMembers;  j-H-) 
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} 


{ . 


} 


evalMemberCrite(TheVar->inember[j],  count,  valid); 
for(k=j-l;  k>=0;  k-) 

compareMemberCrite(TheVar->niember(j], 
TheVar^inember[k],  count,  valid); 

} 


void  evalInMemCrite(C,  count,  valid) 

/*  This  function  will  be  used  to  detennine  the  number  the  of  criteria  used  */ 

/•  to  calculate  the  fitness  ♦/ 

strxict  controller  *C; 

double  ♦count; 

double  ♦valid; 

{ 

inti; 

for(i=0;  i  <  C->numOfInputVar,  1++) 

evaluateCrite(C->inputVars[i],  count,  valid); 

} 

} 


void  evalOutMemCrite(C,  count,  valid) 

/♦  This  function  will  be  used  to  detennine  the  number  the  of  criteria  used  ♦/ 

/♦  to  calculate  the  fitness  ♦/ 
struct  controller  *C; 
double  ♦count; 
double  ♦valid; 

{ 

inti; 

foi(i=0;  i  <  C->numOfOutputVar,  i-H-) 

evaluateCrite(C->outputVars[i],  count,  valid); 

} 

} 

double  evalMemFitness(C) 

/♦  This  function  will  be  us^  to  calculate  the  fitness  of  the  membership  ♦/ 

/♦  fimctions.  The  fitness  is  detennined  by  the  criteria  for  valid  members  ♦/ 
struct  controller  ♦€; 

{ 

double  ruleCount,  validCrite; 
nileCount  =  0.0; 
validCrite  =  0.0; 

evalInMemCrite(C,  &ruleCount,  AvalidCrite); 
evalOutMemCrite(C,  &ruleCount,  &validCrite); 

/♦ 

fprintf(stdout,  "\nThe  number  of  criteria  =  %.ir  ,raleCount); 
fprintf(stdout,  **\nCrite  that  were  valid  =  %.ir.  validCrite); 

♦/ 

return(validCrite/ruleCount); 

} 

void  countMembVert(member,  count) 

/♦  This  function  is  us^  to  count  the  number  of  vertecies  for  a  single  member  ♦/ 
stnict  membeiObject  ♦member, 
int  ♦count; 

{ 

switch(member->type) 
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{ 


case  I1RIANGULAR: 

♦count  =  ♦count  +  3; 
break; 

case  S^FUNCTION: 

♦count  =  ♦count  +  2; 
break; 

case  Z_FUNCTION: 

♦count  =  ♦count  +  2; 
break; 

default 

5)rintf{stderT,"\ncountMembVert:  Error  in  the  case"; 
exit(l); 

} 

} 

void  count VarBit(TheVar,  bitCount) 

/♦  This  function  vdll  be  used  to  determine  the  number  the  of  bits  of  a  */ 

/♦  a  single  variable  ♦/ 
struct  vaiObject  ♦TheVar, 
int  ♦bitCount; 

{ 

intj; 

int  vertCount; 
vertCount  =  0; 

for(j=0;  j  <  TheVar->numOfMembers;  j++) 

{ 

coimtMembVert(TheVar->member|j],  AvertCount); 

} 

♦bitCount  =  •bitCount  +  (veitCount*(TheVar->numBits)); 

} 

int  countBits(C) 

/♦  This  function  will  be  used  to  determine  the  number  of  bits  needed  ♦/ 
struct  controller  ♦C; 

int  bitCount,  i,  j; 
bitCount  =  0; 

for(j=0;j  <  C->numOfInputVar,  j-H-) 

{ 

countVarBit(C->inputVars[jl,  AbitCount); 

} 

for(i=0;  i  <  C^numOfOutputVar,  i-H-) 

{ 

count  VarBit(C->outputVars[i],  AbitCount); 

} 

§}rintf(stdout,  "\nThe  number  of  bits  ==  %d.\n",  bitCount); 
retum(bitCount); 

} 

vertexType  getRealNumber(critter,start,stop4iumBits,lowerBound,delta) 

/♦  This  function  is  used  to  get  real  numbers  from  the  chrome  ♦/ 

unsigned  ♦critter, 

int  ♦start; 

int  ♦stop; 

int  numBits; 

vertexType  lowerBound; 

vertexType  delta; 

{ 

vertexType  realValue; 
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/*  section  of  chromosome  containing  current  integer  field  *! 

♦start  =  ♦stop+  1; 

♦stop  =  ♦start  +  numBits  - 1 ; 

/♦  printf("\n  start  =  %d ,  stop  =  %d  Ichrom  =  %d",  ♦start,  ♦stop, 
Ichrom);  ♦/ 

/♦  check  if  enough  bits  remain,  if  not,  exit  program  ♦/ 
if(^stop  >  Ichrom) 

{ 

§)rintf(stderr,”\nError  in  getRealNumber  fimction.\n"); 
exit(-l); 

} 

realValue  =((double)ithruj2int( ♦start,  ♦stop,  critter))/ 
po\v(2.0,  (doubleXnumBits)); 
realValue  =  reaiValue^delta  +  lowerBound; 
retuni(realValue); 

} 

void  getInValues(critter,  C,  start,  stop) 

/♦  This  function  is  used  to  get  values  for  the  vertecies  of  ♦/ 

/♦  the  input  variables  ♦/ 
unsigned  ♦critter, 
struct  controller  *C\ 
int  ♦start; 
int  ♦stop; 

{ 

intij; 

foT(i=0;  i  <  C->numOfInputVar,  i-H-) 

for(j=0;  j  <  C->inputVars[i].numOfMembers; )++) 

{ 

switch(C->inputVars[i].member[j].type) 

{ 

case  TRIANGULAR: 

C->inputVars[i].memberU].leflVertex  = 
getRealNumber(critter,start,stop,C->inputVars[i],numBits, 
C->inputVars[i].lowerBound,  C->inputVars[i]. delta); 
C->inputVars[i].memb^].centerVertex  = 
getRealNumber(critter,start,stop,C^inputVars[i].numBits, 
C->inputVars[i] .  lowerBound,  C->inputVars[i]  .delta); 
C->inputVars[i].memb^].rightVertex  = 
getRealNumbeifcritter,start,stop,C->inputVars[i].nuinBits, 
C->inputVars[i].lowerBound,  C->inputVars[i]. delta); 

break; 

case  S_FUNCTION: 

C->inputVars[i].memberlj].leftVertex  = 
getReaINumber(critter,start,stop,C->inputVars[i].numBits, 
C->inputVars[i].  lowerBound,  C->inputVars[i]  .delta); 
C->inputVars[i].memb^].centerVertex  = 
getReaINumbeT(critter,start,stop,C->inputVars[i].numBits, 
C^inputVars[i],lowerBound,  C->ii^>utVars[i]. delta); 

break; 

case  Z_FUNCTION: 

C->mputVars[i].member|j].centerVertex  = 
getReaINumber(critter,start,stop,C->mputVars[i].nuinBits, 
C->inputVars[i],lowerBound,  C->inputVars[i].delta); 
C->inputVars[i].memb^].rightVertex  = 
getReaINumb^critter,start,stop,C->inputVars[i].numBits, 
C->inputVars[i]  .lowerBound,  C->inputVars[i].delta); 

break; 

default: 
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void  getValues(critter,  C) 

/♦  This  function  is  used  to  get  the  values  for  the  vertecies.  */ 
unsigned  *critter, 
struct  controller  *0; 

{ 

int  start,  stop; 
start  =  0; 
stop  =  0; 

getInValues(critter,  C,  &start,  &stop); 
getOutValues(critter,  C,  &start,  &stop); 

} 

void  cloneInputVar(original,  clone) 

/*  This  function  is  used  clone  the  input  variables  */ 
struct  controller  *<Higinal; 
struct  controller  *clone; 

{ 

intij; 

for(i=0;  i  <  clone->numOfInputVar,  i-H-) 

{ 

strcpy(clone->input  Vars[i] .  varName, 

original->inputVars[i].  varName); 
clone->inputVars[i]  .numOfMembers= 

original->input  Vars[i]  .numOfMembers; 
clone->inputVars[i].lowerBound= 

original->inputVars[i].lowerBound; 
clone->inputVars[i]  .upperBound= 

original->mputVars[i]  .upperBound; 

clone->inputVars[i].delta= 

original->inputVars[i]  .delta; 
clone->inputVars[i]  .numBits= 

originaI->input  Vars[i]  .numBits; 
if((clone->inputVars[i].member  =  (struct  membeiObject  *) 
malloc((clone->inputVars[i]  .numOfMembers)  * 
sizeof(struct  memberObject)))=NULL) 

{ 

fprintf(stderT,’'not  enough  memory  for  the  clone's "); 
Q)rintf(stderr, "member  functions  of  variable  "); 
fprintf(stderr,"'%s'.\n" ,  clone->inputVars[i] .  varName); 
exit(-l); 

} 

for(j=0;  j  <  clone->inputVars[i]. numOfMembers;  j-H-) 

{ 

strcpy(clone->inputVars[i]  .member{j]  .memName, 

original->ii^)utVars[i].member[j]  .memName); 
clQne->inputVars[i]  jnember[j].type  = 

original->inputVars[i] .  member[j]  .type; 

} 

} 

} 

void  cloneOutputVar(original,  clone) 

/•  This  function  is  used  to  clone  the  output  variables  */ 
struct  controller  *original; 
struct  controller  *clone; 

{ 

inti,j; 

for(i=0;  i  <  clone->numOfOutputVar,  i-H-) 

{ 

strcpy(clone->outputVars[i].varName, 

original->outputVars[i].varName); 
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} 


clone->output  Vars[i]  .numOfMembers= 

original->outputVars[i].numO£Members; 
clone->output  Vars[i] .  lowerBoimd= 

original->output  Vars[i] .  lowerBound; 
clone->output  Vars[i]  .upperBound= 

original->outputVars[i]  .upperBound; 

clone->outputVars[i].delta= 

original->outputVars[i]  .delta; 
clone->outputVars[i].niiinBits= 

origiiial->outputVars[i].numBits; 
if((clone->outputVars[i].member  =  (struct  memberObject  ♦) 
iiialloc((clone->outputVars[i].numOfMembers)  ♦ 
sizeof(struct  memberObject)))=NULL) 

{ 

5)rintf{stderr,’'not  enough  memory  for  the  clone’s "); 
Q)rmtf(stderr,"member  functions  of  variable  ”); 
§>rintf(stderr,'"%s’.\n",  clone->outputVars[i].varName); 
exit(-l); 

} 

far(j=0;j  <clone->outputVars[i].numOfMembers;  j-H-) 

{ 

strcpy(cIone->outputVars[i].member[j],memName, 

original->outputVars[i].member[j].memName); 
clone->outputVars[i].memberU].type  = 

original->outputVar?[i].member[j].type; 


} 


void  cloneController(original,  clone) 

/*  This  function  is  used  to  clone  a  controller  */ 
struct  controller  *original; 
struct  controller  *clone; 

{ 

strcpy(clone->controllerName,  original->controllerName); 
cIone->numOflnputVar  =  original->numOfInputVar, 
clone->numOfOutputVar  =  original->numOfOutputVar, 
if((clone«>inputVars  =  (struct  varObject  *) 

malloc((clone->numOfInputVar)*sizeof(struct  varObject))>=NULL) 

{ 

fprintf(stderr,"not  enough  memory  for  %s.\n", "clone  inputVars"); 
exit(.l); 

} 

if((clone->outputVars  =  (struct  varObject  *) 

malloc((clone->numOfOutputVar)*sizeof(struct  varObject)))==NULL) 

{ 

fprintf(stderr,''not  enough  memory  for  %s.\n’’,"clone  outputVars"); 
exit(-l); 

} 

cloneInputVar(original,  clone); 
cloneOutputVaiforiginal,  clone); 

} 


void  store_data(mt  Generation) 

/♦  This  function  is  used  to  store  the  data  to  weightFile.  The  data  will  be  */ 
/•  used  to  test  the  results  of  the  sga.  */ 

{ 

FILE  *Q>StoreData; 
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int  vertex; 

charWeightFile[120]; 


/♦  File  used  to  store  the  vertex  for  this  ♦/ 

/*  generation.  V 

sprintf(WeightFile,"%s.%d" ,  VertexF  ileBase,Generation); 
if((Q)StoreE)ata=fopen(WeightFile,"w"))  =  NULL) 

{ 

printf("Could  not  open  %s\n”,WeightFile); 
exit(-l); 

} 

fprintf(fpStoreData,"%s\n"JDataFile); 

§irintf(^toreData,"%d\n"J4umbeiOfSamples); 

^printf(j^StoreData,"%lf\n'’,  bestfitfitness); 

Q)rintf(fi)StoreData,  "%d\n'', Generation); 

Q)rintf(fi)StoreE>ata,  **%d\n"4iiaxgen); 

Q>rintf(Q>StoreData,"%f\n**,  pcross); 
fprintf(f^toreData,’'%f\n",  pmutation); 

5MTntf(^toreData,  '’%d\n",popsize); 
fclose(^toreData); 

}/*  End  of  function  store^data.  ♦/ 


application() 

/*  This  routine  should  contain  any  application-dependent  */ 

/•  computations  that  should  be  perfonned  before  eart'  GA  cycle  ♦/ 
/•  called  by  main()  */ 

{ 

} 


app_data() 

/*  application  dependent  data  input,  called  by  init_data()  ♦/ 

/*  ask  your  input  questions  here,  and  put  output  in  global  variables  •/ 

{ 

} 


app_free() 

/*  This  routine  should  free  any  memory  allocated  */ 

/*  in  the  appbcation-dependent  routines,  called  by  freeall()  */ 

{ 

intj; 

freeVariables(&person);  /•  deallocate  variable  space  */ 
freeRules(&ruleBase);  /•  deallocate  rule  space  •/ 

/*  deallocate  space  for  data  variables  */ 
for(j  =  0;  j  <  NumberOfSamples;  j-H-) 

{ 

free(InputData[j]); 

free(OutputData[j]); 

fiee(CalculatedOutput[j]); 

} 

free(InputData); 

fi:ee(C)u^utI>ata); 

fi:iee(CalculatedOutput); 

} 

app_init() 

/*  Application  dependent  initialization  routine  called  by  initialize().  */ 

{ 

} 
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APPENDIX  D 

SAMPLE  RUN  OF  THE  GA  TUNER 


I  SGA-C  (vl .  1)  -  A  Simple  Genetic  Algorithm  | 

I  (c)  David  E.  Goldberg  1986,  All  Rights  Reserved  | 

I  C  version  by  Robert  E.  Smith,  U.  of  Alabama  | 

I  vl.l  modifications!^  Jeff Earickson,  Boeing  Company  | 


SGA  Parameters 


Total  Population  size  =  50 

Chromosome  length  (Ichrom)  =  210 

Maximum  #  of  generations  (maxgen)  =  200 

Crossover  probability  (paoss)  =  0.600000 
Mutation  probability  (pmutation)  =  0.010000 


RUN  1  of  1;  GENERATION  4->200 


Generation  4  Accumulated  Statistics: 

Total  Crossovers  =  80,  Total  Mutations  =  529 

min  =  0.347222  max  =  0.430556  avg  =  0.402778  sum  =  20.138889 

Gl(^  Best  Individual  so  far.  Generation  4: 

Fitness  =  0.430556: 


CONTROLLER:  simplePend 
- INPUTS - 


number  of  input  varibles  =  2 
variable[l]  =  'error*  and  it  has  3  members: 


Z-Function:  member[l]  =  negative 

centerVertex  =  -0.850 
rightVertex  =0.711 

Tiangular:  member[2]  =  zero 

leftVertex  =  -0.938 
centerVertex  =  0. 195 
rightVertex  =0.771 

S-Function:  member[3]  =  positive 

leftVertex  =0.557 
centerVertex  =  0.71 1 
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lower  bounds  =  -1.000000 
upper  bounds  =  1.000000 
delta  =2.000000 
num.  of  bits  =  10 

variable[2]  =  Derror'  and  it  has  3  members: 

Z-Function;  meinber[l]  =  negative 

centerVertex  =  0.527 
rightVertex  =-0.809 

Tiangular:  member[2]  =  zero 

IcftVertex  =-0.498 
centerVertex  =  -0.496 
rightVertex  =0.090 

S-Function;  member[3]  =  positive 

leftVertex  =  -0.029 
centerVertex  =  0.947 


lower  bounds  =  -1.000000 
upper  bounds  =  1.000000 
delta  =  2.000000 
num.  of  bits  =  10 


- OUTPUTS - 

number  of  output  varibles  =  1 

variable[l]  =  'current'  and  it  has  3  members: 

Z-Function:  member[l]  =  negative 

centerVertex  =  -0.240 
rightVertex  =0.586 

Tiangular:  member[2]  =  zero 

leftVertex  =  -0.967 
centerVertex  =  -0.848 
rightVertex  =0.029 

S-Function:  member[3]  =  positive 

leftVertex  =-0.908 
centerVertex  =  -0.082 


lower  bounds  =  -1.000000 
upper  bounds  =  1 .000000 
delta  =  2.000000 
num.  of  bits  =  10 


RUN  1  of  1:  GENERATION  29->200 

Generation  29  Accumulated  Statistics; 

Total  Crossovers  =  464,  Total  Mutations  =  3085 
min  =  0.458333  max  =  0.9623 14  avg  =  0.795700 
Global  Best  Individual  so  far.  Generation  28: 

Fitness  =  0.963756: 

sum  =  39.784985 

RUN  1  of  1:  GENERATION  59->200 

k  i-* 

Generation  59  Accumulated  Statistics: 

Total  Crossovers  =  903,  Total  Mutations  =  6188 
min  =  0.430556  max  =  0.985419  avg  =  0.808049 
Global  Best  Individual  so  far.  Generation  59: 

Fitness  =  0.985419: 

sum  =  40.402465 

RUN  1  of  1:  GENERATION  79->200 

Generation  79  Accuniulated  Statistics: 

Total  Crossovers  =  1205,  Total  Mutations  =  8258 
min  -  0.444444  max  =  0.989729  avg  =  0.826472 
Global  Best  Individual  so  far.  Generation  78: 

Fitness  =  0.989907: 

sum  =  41.323584 

RUN  1  of  1:  GENERATION  99->200 

Generation  99  Accumulated  Statistics: 

Total  Crossovers  =  1509,  Total  Mutations  =  10320 
min  =  0.430556  max  =  0.991671  avg  =  0.776286 
Global  Best  Individual  so  far.  Generation  97: 

Fitness  =  0.991671: 

sum  =  38.814297 

CONTROLLER:  simplePend 

- INPUTS - 

number  of  input  varibles  =  2 
variable[l]  =  'error'  and  it  has  3  members: 

Z-Function:  member[l]  =  negative 
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centerVertex  =  -0.650 
rightVertex  =0.043 

Tiangular:  inember[2]  =  zero 

leftVertex  =-0.389 
centerVertex  =  0.025 
rightVertex  =0.502 

S-Function:  member[3]  =  positive 

leftVertex  =-0.012 
centerVertex  =  0.3'^9 


lower  bounds  =  -1.000000 
upper  bounds  =  1 .000000 
delta  =2.000000 
num.  of  bits  =  10 

variable[2]  =  Derror*  and  it  has  3  members: 

Z-Function:  member[l]  =  negative 

centerVertex  =  -0.273 
rightVertex  =0.018 

Tiangular:  member[2]  =  zero 

leftVertex  =-0.451 
centerVertex  =  -0. 199 
rightVertex  =0.502 

S-Function:  member[3]  =  positive 

leftVertex  =  -0.035 
centerVertex  =  0.562 


lower  bounds  =  -1.000000 
upper  bounds  =  1.000000 
delta  =  2.000000 
num.  of  bits  =  10 


- OUTPUTS - 

number  of  output  varibles  =  1 
variable[l]  =  'current'  and  it  has  3  members: 

Z-Function:  member[l]  =  negative 
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centerVertex  =  -0.61 1 
rightVertex  =0.047 


Tiangular;  inember[2]  =  zero 

left  Vertex  =-0.297 
centerVertex  =  -0. 178 
rightVertex  =0.471 

S-Function:  member[3]  =  positive 

left  Vertex  =0.102 
centerVertex  =  0.369 


lower  bounds  =  -1.000000 
upper  bounds  =  1.000000 
delta  =2.000000 
num.  of  bits  =  10 


RUN  1  of  1:  GENERATION  199->200 


Generation  199  Accumulated  Statistics; 

Total  Crossovers  =  3000,  Total  Mutations  =  20714 

min  =  0  444444  max  =  0.992721  avg  =  0.812836  sum  =  40.641800 

Global  Best  Individual  so  far.  Generation  191: 

Fitness  =  0.992738; 
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Appendix  E 

Rule  Space  for  Pendulum  Balancing  Controller 

IF  error  IS  negative  THEN  current  IS  negative  END 
IF  error  IS  negative  THEN  current  IS  zero  END 
IF  error  IS  negative  THEN  current  IS  positive  END 

IF  error  IS  negative  AND  Derror  IS  negative  THEN  current  IS  negative  END 
IF  error  IS  negative  AND  Derror  IS  negative  THEN  current  IS  zero  END 
IF  error  IS  negative  AND  Derror  IS  negative  THEN  current  IS  positive  END 
IF  error  IS  negative  AND  Derror  IS  zero  THEN  current  IS  negative  END 
IF  error  IS  negative  AND  Derror  IS  zero  THEN  current  IS  zero  END 
IF  error  IS  negative  AND  Derror  IS  zero  THEN  current  IS  positive  END 
IF  error  IS  negative  AND  Derror  IS  positive  THEN  current  IS  negative  END 
IF  error  IS  negative  AND  Derror  IS  positive  THEN  current  IS  zero  END 
IF  error  IS  negative  AND  Derror  IS  positive  THEN  current  IS  positive  END 
IF  error  IS  zero  THEN  current  IS  negative  END 
IF  error  IS  zero  THEN  current  IS  zero  END 
IF  error  IS  zero  THEN  current  IS  positive  END 

IF  error  IS  zero  AND  Derror  IS  negative  THEN  current  IS  negative  END 
IF  error  IS  zero  AND  Derror  IS  negative  THEN  current  IS  zero  END 
IF  error  IS  zero  AND  Derror  IS  negative  THEN  current  IS  positive  END 
IF  error  IS  zero  AND  Derror  IS  zero  THEN  current  IS  negative  END 
IF  error  IS  zero  AND  Derror  IS  zero  THEN  current  IS  zero  END 
IF  error  IS  zero  AND  Derror  IS  zero  THEN  current  IS  positive  END 
IF  error  IS  zero  AND  Derror  IS  positive  THEN  current  IS  negative  END 
IF  error  IS  zero  AND  Derror  IS  positive  THEN  current  IS  zero  END 
IF  error  IS  zero  AND  Derror  IS  positive  THEN  current  IS  positive  END 
IF  error  IS  positive  THEN  current  IS  negative  END 
IF  error  IS  positive  THEN  current  IS  zero  END 
IF  error  IS  positive  THEN  current  IS  positive  END 

IF  error  IS  positive  AND  Derror  IS  negative  THEN  current  IS  negative  END 
IF  error  IS  positive  AND  Derror  IS  negative  THEN  current  IS  zero  END 
IF  error  IS  positive  AND  Derror  IS  negative  THEN  current  IS  positive  END 
IF  error  IS  positive  AND  Derror  IS  zero  THEN  current  IS  negative  END 
IF  error  IS  positive  AND  Derror  IS  zero  THEN  current  IS  zero  END 
IF  error  IS  positive  AND  Derror  IS  zero  THEN  current  IS  positive  END 
IF  error  IS  positive  AND  Derror  IS  positive  THEN  current  IS  negative  END 
IF  error  IS  positive  AND  Derror  IS  positive  THEN  current  IS  zero  END 
IF  error  IS  positive  AND  Derror  IS  positive  THEN  current  IS  positive  END 
IF  Derror  IS  negative  THEN  current  IS  negative  END 
IF  Derror  IS  negative  THEN  current  IS  zero  END 
IF  Derror  IS  negative  THEN  current  IS  positive  END 
IF  Derror  IS  zero  THEN  current  IS  negative  END 
IF  Derror  IS  zero  THEN  current  IS  zero  END 
IF  Derror  IS  zero  THEN  current  IS  positive  END 
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Derror  IS  positive  THEN  current  IS  negative  END 
Derror  IS  positive  THEN  current  IS  zero  END 
Derror  IS  positive  THEN  current  IS  positive  END 
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Appendix  F 

Rules  Evolved  by  the  Genetic  Algorithm 
IF  error  IS  negative  THEN  current  IS  positive  END 
IF  error  IS  zero  AND  Derror  IS  negative  THEN  current  IS  positive  END 
IF  error  IS  zero  AND  Derror  IS  zero  THEN  current  IS  zero  END 
IF  error  IS  positive  THEN  current  IS  negative  END 
IF  error  IS  positive  AND  Derror  IS  negative  THEN  current  IS  positive  END 
IF  Derror  IS  positive  THEN  current  IS  negative  END 


CONTROLLER:  simplePenduIum 

- INPUTS - - 

number  of  input  varibles  =  2 

variable[l]  =  ’error’  and  it  has  3  members: 

Z-Function:  member[l]  =  negative 

center  Vertex  =  -0.443 
rightVertex  =  -0. 1 72 

Tiangular:  member[2]  =  zero 

leftVertex  =  -0.955 
centerVertex  =  0.291 
rightVertex  =0.301 

S-Function:  member[3]  =  positive 

leftVertex  =-0.014 
centerVertex  =  0.938 


lower  bounds  =  -1.000000 
upper  bounds  =  1.000000 
delta  =  2.000000 
num.  of  bits  =  10 
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variable[2]  =  Derror’  and  it  has  3  members: 

Z-Function:  member[l]  =  negative 

centerVertex  =  -0.504 
rightVertex  =-0.172 

Tiangular:  member[2]  =  zero 

leftVertex  =  -0.365 
centerVertex  =  -0.244 
rightVertex  =0.746 

S-Function:  member[3]  =  positive 

leftVertex  =-0.014 
centerVertex  =  0.680 


lower  bounds  =  -1.000000 
upper  bounds  =  1 .000000 
delta  =  2.000000 
num.  of  bits  =  10 


- OUTPUTS - 

number  of  output  varibles  =  1 

variable!  1]  =  ’current’  and  it  has  3  members: 

Z-Function:  member!  1  ]  =  negative 

centerVertex  =  -0.283 
rightVertex  =-0.152 

Tiangular:  member!2]  =  zero 

leftVertex  =-0.016 
centerVertex  =  0. 102 
rightVertex  =0.256 

S-Function:  member!3]  =  positive 
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leftVertex  =-0.012 
centerVertex  =  0.566 


lower  bounds  =  -1.000000 
upper  bounds  =  1.000000 
delta  =  2.000000 
num.  of  bits  =  10 


The  number  of  criteria  =  36 
Crite  that  were  valid  =  36 

Mean  Square  Error  =  0.039 
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Appendix  G 

Rule  Space  for  Tank  Recognition  Problem 

ini  IS  low  THEN  class  IS  low  END 

ini  IS  low  THEN  class  IS  med  END 

ini  IS  low  THEN  class  IS  high  END 

ini  IS  low  AND  in2  IS  low  THEN  class  IS  low  END 

in  1  IS  low  AND  in2  IS  low  THEN  class  IS  med  END 

in  1  IS  low  AND  in2  IS  low  THEN  class  IS  high  END 

in  1  IS  low  AND  in2  IS  low  AND  in3  IS  low  THEN  class  IS  low  END 

ini  IS  low  AND  in2  IS  low  AND  in3  IS  low  THEN  class  IS  med  END 

in  1  IS  low  AND  in2  IS  low  AND  in3  IS  low  THEN  class  IS  high  END 

in  I  IS  low  AND  in2  IS  low  AND  in3  IS  med  THEN  class  IS  low  END 

ini  IS  low  AND  in2  IS  low  AND  in3  IS  med  THEN  class  IS  med  END 

ini  IS  low  AND  in2  IS  low  AND  in3  IS  med  THEN  class  IS  high  END 

ini  IS  low  AND  in2  IS  low  AND  in3  IS  high  THEN  class  IS  low  END 

ini  IS  low  AND  in2  IS  low  AND  in3  IS  high  THEN  class  ISmed  END 

ini  IS  low  AND  in2  IS  low  AND  in3  IS  high  THEN  class  IS  high  END 

ini  IS  low  AND  in2  IS  med  THEN  class  IS  low  END 

ini  IS  low  AND  in2  IS  med  THEN  class  IS  med  END 

ini  IS  low  AND  in2  IS  med  THEN  class  IS  high  END 

ini  IS  low  AND  in2  IS  med  AND  in3  IS  low  THEN  class  IS  low  END 

in  1  IS  low  AND  in2  IS  med  AND  in3  IS  low  THEN  class  IS  med  END 

inl  IS  low  AND  in2  IS  med  AND  in3  IS  low  THEN  class  IS  high  END 

nl  IS  low  AND  in2  IS  med  AND  in3  IS  med  THEN  class  IS  low  END 

n  1  IS  low  AND  in2  IS  med  AND  in3  IS  med  THEN  class  IS  med  END 

nl  IS  low  AND  in2  IS  med  AND  in3  IS  med  THEN  class  IS  high  END 

n  1  IS  low  AND  in2  IS  med  AND  in3  IS  high  THEN  class  IS  low  END 

nl  IS  low  AND  in2  IS  med  AND  in3  IS  high  THEN  class  ISmed  END 

n  1  IS  low  AND  in2  IS  med  AND  in3  IS  high  THEN  class  IS  high  END 

nl  IS  low  AND  in2  IS  high  THEN  class  IS  low  END 

nl  IS  low  AND  in2  IS  high  THEN  class  ISmed  END 

nl  IS  low  AND  in2  IS  high  THEN  class  IS  high  END 

nl  IS  low  AND  in2  IS  high  AND  in3  IS  low  THEN  class  IS  low  END 

nl  IS  low  AND  in2  IS  high  AND  in3  IS  low  THEN  class  ISmed  END 

nl  IS  low  AND  in2  IS  high  AND  in3  IS  low  THEN  class  IS  high  END 

nl  IS  low  AND  in2  IS  high  AND  in3  ISmed  THEN  class  IS  low  END 

nl  IS  low  AND  in2  IS  high  AND  in3  ISmed  THEN  class  IS  med  END 

nl  IS  low  AND  in2  IS  high  AND  in3  ISmed  THEN  class  IS  high  END 

nl  IS  low  AND  in2  IS  high  AND  in3  IS  high  THEN  class  IS  low  END 

nl  IS  low  AND  in2  IS  high  AND  in3  IS  high  THEN  class  ISmed  END 

nl  IS  low  AND  in2  IS  high  AND  in3  IS  high  THEN  class  IS  high  END 

nl  IS  low  AND  in3  IS  low  THEN  class  IS  low  END 

nl  IS  low  AND  in3  IS  low  THEN  class  IS  med  END 

nl  IS  low  AND  in3  IS  low  THEN  class  IS  high  END 
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IF  ini  IS  low  AND  in3  IS  med  THEN  class  IS  low  END 

EF  ini  IS  low  AND  in3  IS  med  THEN  class  IS  med  END 

IF  ini  IS  low  AND  in3  IS  med  THEN  class  IS  high  END 

IF  ini  IS  low  AND  in3  IS  high  THEN  class  IS  low  END 

IF  ini  IS  low  AND  in3  IS  high  THEN  class  IS  med  END 

IF  ini  IS  low  AND  in3  IS  high  THEN  class  IS  high  END 

IF  in  I  IS  med  THEN  class  IS  low  END 

IF  ini  IS  med  THEN  class  IS  med  END 

IF  in  1  IS  med  THEN  class  IS  high  END 

IF  in  1  IS  med  AND  in2  IS  low  THEN  class  IS  low  END 

IF  ini  IS  med  AND  in2  IS  low  THEN  class  IS  med  END 

IF  ini  IS  med  AND  in2  IS  low  THEN  class  IS  high  END 

IF  in  I  IS  med  AND  in2  IS  low  AND  in3  IS  low  THEN  class  IS  low  END 

IF  ini  IS  med  AND  in2  IS  low  AND  in3  IS  low  THEN  class  IS  med  END 

IF  ini  IS  med  AND  in2  IS  low  AND  in3  IS  low  THEN  class  IS  high  END 

IF  ini  IS  med  AND  in2  IS  low  AND  in3  IS  med  THEN  class  IS  low  END 

IF  ini  IS  med  AND  in2  IS  low  AND  in3  IS  med  THEN  class  IS  med  END 

IF  ini  IS  med  AND  in2  IS  low  AND  in3  IS  med  THEN  class  IS  high  END 

IF  ini  IS  med  AND  in2  IS  low  AND  in3  IS  high  THEN  class  IS  low  END 

IF  ini  IS  med  AND  in2  IS  low  AND  in3  high  THEN  class  IS  med  END 

IF  ini  IS  med  AND  in2  IS  low  AND  in3  IS  high  THEN  class  IS  high  END 

IF  in  1  IS  med  AND  in2  IS  med  THEN  class  IS  low  END 

IF  in  1  IS  med  AND  in2  IS  med  THEN  class  IS  med  END 

IF  ini  IS  med  AND  in2  IS  med  THEN  class  IS  high  END 

IF  ini  IS  med  AND  in2  IS  med  AND  in3  IS  low  THEN  class  IS  low  END 

IF  ini  IS  med  AND  in2  IS  med  AND  in3  IS  low  THEN  class  IS  med  END 

IF  ini  IS  med  AND  in2  IS  med  AND  in3  IS  low  THEN  class  IS  high  END 

IF  ini  IS  med  AND  in2  IS  med  AND  in3  IS  med  THEN  class  IS  low  END 

IF  ini  IS  med  AND  in2  IS  med  AND  in3  IS  med  THEN  class  IS  med  END 

EF  inl  IS  med  AND  in2  IS  med  AND  in3  IS  med  THEN  class  IS  high  END 

IF  ini  IS  med  AND  in2  IS  med  AND  in3  IS  high  THEN  class  IS  low  END 

IF  inl  IS  med  AND  in2  IS  med  AND  in3  IS  high  THEN  class  IS  med  END 

IF  inl  IS  med  AND  in2  IS  med  AND  in3  IS  high  THEN  class  IS  high  END 

IF  inl  IS  med  AND  in2  IS  high  THEN  class  IS  low  END 

IF  inl  IS  med  AND  in2  IS  high  THEN  class  IS  med  END 

IF  inl  IS  med  AND  in2  IS  high  THEN  class  IS  high  END 

EF  inl  IS  med  AND  in2  IS  high  AND  in3  IS  low  THEN  class  IS  low  END 

IF  inl  IS  med  AND  in2  IS  high  AND  in3  IS  low  THEN  class  ISmed  END 

IF  inl  IS  med  AND  in2  IS  high  AND  in3  IS  low  THEN  class  IS  high  END 

IF  inl  IS  med  AND  in2  IS  high  AND  in3  ISmed  THEN  class  IS  low  END 

IF  inl  IS  med  AND  in2  IS  high  AND  in3  ISmed  THEN  class  IS  med  END 

IF  inl  IS  med  AND  in2  IS  high  AND  in3  ISmed  THEN  class  IS  high  END 

IF  inl  IS  med  AND  in2  IS  high  AND  in3  IS  high  THEN  class  IS  low  END 

IF  inl  IS  med  AND  in2  IS  high  AND  in3  IS  high  THEN  class  ISmed  END 

IF  inl  IS  med  AND  in2  IS  high  AND  in3  IS  high  THEN  class  IS  high  END 
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ini  IS  med  AND  in3  IS  low  THEN  class  IS  low  END 
ini  IS  med  AND  in3  IS  low  THEN  class  IS  med  END 
ini  IS  med  AND  in3  IS  low  THEN  class  IS  high  END 
in  1  IS  med  AND  in3  IS  med  THEN  class  IS  low  END 
in  1  IS  med  AND  in3  IS  med  THEN  class  IS  med  END 


ini  IS  med  AND  in3  IS  med  THEN  class  IS  high  END 
ini  IS  med  AND  in3  IS  high  THEN  class  IS  low  END 
ini  IS  med  AND  in3  IS  high  THEN  class  IS  med  END 
ini  IS  med  AND  in3  IS  high  THEN  class  IS  high  END 
ini  IS  high  THEN  class  IS  low  END 
in  1  IS  high  THEN  class  IS  med  END 


il  IS  high  THEN  class  IS  high  END 
il  IS  high  AND  in2  IS  low  THEN  class  IS  low  END 
il  IS  high  AND  in2  IS  low  THEN  class  IS  med  END 
il  IS  high  AND  in2  IS  low  THEN  class  IS  high  END 
Il  IS  high  AND  in2  IS  low  AND  in3  IS  low  THEN  class  IS  low  END 
il  IS  high  AND  in2  IS  low  AND  in3  IS  low  THEN  class  ISmed  END 
il  IS  high  AND  in2  IS  low  AND  in3  IS  low  THEN  class  IS  high  END 
il  IS  high  AND  in2  IS  low  AND  in3  ISmed  THEN  class  IS  low  END 
il  IS  high  AND  in2  IS  low  AND  in3  ISmed  THEN  class  IS  med  END 
il  IS  high  AND  in2  IS  low  AND  in3  ISmed  THEN  class  IS  high  END 
1  IS  high  AND  in2  IS  low  AND  in3  IS  high  THEN  class  IS  low  END 
1  IS  high  AND  in2  IS  low  AND  in3  IS  high  THEN  class  ISmed  END 
1  IS  high  AND  in2  IS  low  AND  in3  IS  high  THEN  class  IS  high  END 
1  IS  high  AND  in2  ISmed  THEN  class  IS  low  END 
1  IS  high  AND  in2  IS  med  THEN  class  IS  med  END 
1  IS  high  AND  in2  ISmed  THEN  class  IS  high  END 


1  IS  high  AND  in2  ISmed  AND  in3  IS  low  THEN  class  IS  low  END 
1  IS  high  AND  in2  ISmed  AND  in3  IS  low  THEN  class  IS  med  END 
1  IS  high  AND  in2  ISmed  AND  in3  IS  low  THEN  class  IS  high  END 
1  IS  high  AND  in2  ISmed  AND  in3  IS  med  THEN  class  IS  low  END 
1  IS  high  AND  in2  ISmed  AND  in3  IS  med  THEN  class  IS  med  END 
1  IS  high  AND  in2  ISmed  AND  in3  IS  med  THEN  class  IS  high  END 
1  IS  high  AND  in2  ISmed  AND  in3  IS  high  THEN  class  IS  low  END 
1  IS  high  AND  in2  ISmed  AND  in3  IS  high  THEN  class  ISmed  END 
1  IS  high  AND  in2  ISmed  AND  in3  IS  high  THEN  class  IS  high  END 
1  IS  high  AND  in2  IS  high  THEN  class  IS  low  END 
1  IS  high  AND  in2  IS  high  THEN  class  ISmed  END 
1  IS  high  AND  in2  IS  high  THEN  class  IS  high  END 
1  IS  high  AND  in2  IS  high  AND  in3  IS  low  THEN  class  IS  low  END 
1  IS  high  AND  in2  IS  high  AND  in3  IS  low  THEN  class  ISmed  END 
1  IS  high  AND  in2  IS  high  AND  in3  IS  low  THEN  class  IS  high  END 
1  IS  high  AND  in2  IS  high  AND  in3  ISmed  THEN  class  IS  low  END 
1  IS  high  AND  in2  IS  high  AND  in3  ISmed  THEN  class  IS  med  END 
1  IS  high  AND  in2  IS  high  AND  in3  ISmed  THEN  class  IS  high  END 
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IF  ini  IS  high  AND  in2  IS  high  AND  in3  IS  high  THEN  class  IS  low  END 

IF  ini  IS  high  AND  in2  IS  high  AND  in3  IS  high  THEN  class  ISmed  END 

IF  ini  IS  high  AND  in2  IS  high  AND  in3  IS  high  THEN  class  IS  high  END 

IF  ini  IS  high  AND  in3  IS  low  THEN  class  IS  low  END 

IF  ini  IS  high  AND  in3  IS  low  THEN  class  ISmed  END 

IF  ini  IS  high  AND  in3  IS  low  THEN  class  IS  high  END 

IF  ini  IS  high  AND  in3  ISmed  THEN  class  IS  low  END 

IF  ini  IS  high  AND  in3  ISmed  THEN  class  IS  med  END 

IF  in  I  IS  high  AND  in3  ISmed  THEN  class  IS  high  END 

IF  ini  IS  high  AND  in3  IS  high  THEN  class  IS  low  END 

IF  ini  IS  high  AND  in3  IS  high  THEN  class  ISmed  END 

IF  ini  IS  high  AND  in3  IS  high  THEN  class  IS  high  END 

IF  in2  IS  low  THEN  class  IS  low  END 

IF  in2  IS  low  THEN  class  IS  med  END 

IF  in2  IS  low  THEN  class  IS  high  END 

IF  in2  IS  low  AND  in3  IS  low  THEN  class  IS  low  END 

IF  in2  IS  low  AND  in3  IS  low  THEN  class  IS  med  END 

IF  in2  IS  low  AND  in3  IS  low  THEN  class  IS  high  END 

IF  in2  IS  low  AND  in3  IS  med  THEN  class  IS  low  END 

IF  in2  IS  low  AND  in3  IS  med  THEN  clas^  IS  med  END 

IF  in2  IS  low  AND  in3  IS  med  THEN  class  IS  high  END 

IF  in2  IS  low  AND  in3  IS  high  THEN  class  IS  low  END 

IF  in2  IS  low  AND  in3  IS  high  THEN  class  IS  med  END 

IF  in2  IS  low  AND  in3  IS  high  THEN  class  IS  high  END 

IF  in2  IS  med  THEN  class  IS  low  END 

IF  in2  IS  med  THEN  class  IS  med  END 

IF  in2  IS  med  THEN  class  IS  high  END 

IF  in2  IS  med  AND  in3  IS  low  THEN  class  IS  low  END 

IF  in2  IS  med  AND  in3  IS  low  THEN  class  IS  med  END 

IF  in2  IS  med  AND  in3  IS  low  THEN  class  IS  high  END 

IF  in2  IS  med  AND  in3  IS  med  THEN  class  IS  low  END 

IF  in2  IS  med  AND  in3  IS  med  THEN  class  IS  med  END 

IF  in2  IS  med  AND  in3  IS  med  THEN  class  IS  high  END 

IF  in2  IS  med  AND  in3  IS  high  THEN  class  IS  low  END 

IF  in2  IS  med  AND  in3  IS  high  THEN  class  ISmed  END 

IF  in2  IS  med  AND  in3  IS  high  THEN  class  IS  high  END 

IF  in2  IS  high  THEN  class  IS  low  END 

IF  in2  IS  high  THEN  class  IS  med  END 

IF  in2  IS  high  THEN  class  IS  high  END 

IF  in2  IS  high  AND  in3  IS  low  THEN  class  IS  low  END 

IF  in2  IS  high  AND  in3  IS  low  THEN  class  ISmed  END 

IF  in2  IS  high  AND  in3  IS  low  THEN  class  IS  high  END 

IF  in2  IS  high  AND  in3  ISmed  THEN  class  IS  low  END 

IF  in2  IS  high  AND  in3  ISmed  THEN  class  IS  med  END 

IF  in2  IS  high  AND  in3  ISmed  THEN  class  IS  high  END 
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in2  IS  , high  AND  in3  IS  high  THEN  class  IS  low  END 

in2  IS  high  AND  in3  IS  high  THEN  class  ISmed  END 

in2  IS  high  AND  in3  IS  high  THEN  class  IS  high  END 

in3  IS  low  THEN  class  IS  low  END 

in3  IS  low  THEN  class  IS  med  END 

in3  IS  low  THEN  class  IS  high  END 

in3  IS  med  THEN  class  IS  low  END 

in3  IS  med  THEN  class  IS  med  END 

in3  IS  med  THEN  class  IS  high  END 

in3  IS  high  THEN  class  IS  low  END 

in3  IS  high  THEN  class  ISmed  END 

in3  IS  high  THEN  class  IS  high  END 
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Appendix  H 

Rules  Evolved  by  the  Genetic  Algorithm 


IF  ini  IS  negative  AND  in3  IS  negative  THEN  class  IS  positive  END 
IF  ini  IS  positive  THEN  class  IS  positive  END 
IF  in2  IS  negative  AND  in3  IS  negative  THEN  class  IS  negative  END 
IF  in2  IS  negative  AND  in3  IS  negative  THEN  class  IS  positive  END 
IF  in2  IS  positive  AND  in3  IS  negative  THEN  class  IS  negative  END 
IF  in3  IS  negative  THEN  class  IS  negative  END 
IF  in3  IS  positive  THEN  class  IS  positive  END 

Mean  Square  Error  =  0.414 

Floor-Ceiling  Mean  Square  Error  =  0.498 

Floor-Ceiling  Percent  Correctly  Classified  =  394/516  =  76.36 
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